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Foreword 



More than a generation of German-speaking students around the world have worked their 
way to an understanding and appreciation of the power and beauty of modern theoretical 
physics — with mathematics, the most fundamental of sciences — using Walter Greiner’s 
textbooks as their guide. 

The idea of developing a coherent, complete presentation of an entire field of science in a 
series of closely related textbooks is not a new one. Many older physicians remember with 
real pleasure their sense of adventure and discovery as they worked their ways through the 
classic series by Sommerfeld, by Planck, and by Landau and Lifshitz. From the students’ 
viewpoint, there are a great many obvious advantages to be gained through the use of 
consistent notation, logical ordering of topics, and coherence of presentation; beyond this, 
the complete coverage of the science provides a unique opportunity for the author to convey 
his personal enthusiasm and love for his subject. 

These volumes on classical physics, finally available in English, complement Greiner’s 
texts on quantum physics, most of which have been available to English-speaking audiences 
for some time. The complete set of books will thus provide a coherent view of physics that 
includes, in classical physics, thermodynamics and statistical mechanics, classical dynam- 
ics, electromagnetism, and general relativity; and in quantum physics, quantum mechanics, 
symmetries, relativistic quantum mechanics, quantum electro- and chromodynamics, and 
the gauge theory of weak interactions. 

What makes Greiner’s volumes of particular value to the student and professor alike is 
their completeness. Greiner avoids the all too common “it follows that . . . ,” which conceals 
several pages of mathematical manipulation and confounds the student. He does not hesitate 
to include experimental data to illuminate or illustrate a theoretical point, and these data, 
like the theoretical content, have been kept up to date and topical through frequent revision 
and expansion of the lecture notes upon which these volumes are based. 

Moreover, Greiner greatly increases the value of his presentation by including something 
like one hundred completely worked examples in each volume. Nothing is of greater 
importance to the student than seeing, in detail, how the theoretical concepts and tools 



V 




FOREWORD 



under study are applied to actual problems of interest to working physicists. And, finally, 
Greiner adds brief biographical sketches to each chapter covering the people responsible 
for the development of the theoretical ideas and/or the experimental data presented. It 
was Auguste Comte (1789-1857) in his Positive Philosophy who noted, “To understand a 
science it is necessary to know its history.” This is all too often forgotten in modern physics 
teaching, and the bridges that Greiner builds to the pioneering figures of our science upon 
whose work we build are welcome ones. 

Greiner’s lectures, which underlie these volumes, are internationally noted for their 
clarity, for their completeness, and for the effort that he has devoted to making physics an 
integral whole. His enthusiasm for his sciences is contagious and shines through almost 
every page. 

These volumes represent only a part of a unique and Herculean effort to make all of 
theoretical physics accessible to the interested student. Beyond that, they are of enormous 
value to the professional physicist and to all others working with quantum phenomena. 
Again and again, the reader will find that, after dipping into a particular volume to review a 
specific topic, he or she will end up browsing, caught up by often fascinating new insights 
and developments with which he or she had not previously been familiar. 

Having used a number of Greiner’s volumes in their original German in my teaching 
and research at Yale, I welcome these new and revised English translations and would 
recommend them enthusiastically to anyone searching for a coherent overview of physics. 

D. Allan Bromley 

Henry Ford II Professor of Physics 

Yale University 

New Haven, Connecticut, USA 




Preface 



Theoretical physics has become a many faceted science. For the young student, it is difficult 
enough to cope with the overwhelming amount of new material that has to be learned, 
let alone obtain an overview of the entire field, which ranges from mechanics through 
electrodynamics, quantum mechanics, field theory, nuclear and heavy-ion science, statistical 
mechanics, thermodynamics, and solid-state theory to elementary-particle physics; and this 
knowledge should be acquired in just eight to ten semesters, during which, in addition, a 
diploma or master’s thesis has to be worked on or examinations prepared for. All this can be 
achieved only if the university teachers help to introduce the student to the new disciplines 
as early as possible, in order to create interest and excitement that in turn set free essential 
new energy. 

At the Johann Wolfgang Goethe University in Frankfurt am Main, we therefore con- 
front the student with theoretical physics immediately, in the first semester. Theoretical 
Mechanics I and II, Electrodynamics, and Quantum Mechanics I — An Introduction are the 
courses during the first two years. These lectures are supplemented with many mathemati- 
cal explanations and much support material. After the fourth semester of studies, graduate 
work begins, and Quantum Mechanics II — Symmetries, Statistical Mechanics and Ther- 
modynamics, Relativistic Quantum Mechanics, Quantum Electrodynamics, Gauge Theory 
of Weak Interactions, and Quantum Chromodynamics are obligatory. Apart from these, 
a number of supplementary courses on special topics are offered, such as Hydrodynam- 
ics, Classical Field Theory, Special and General Relativity, Many-Body Theories, Nuclear 
Models, Models of Elementary Particles, and Solid-State Theory. 

This volume of lectures. Classical Mechanics: Point Particles and Relativity , deals with 
the first and more elementary part of the important field of classical mechanics. We have 
tried to present the subject in a manner that is both interesting to the student and easily 
accessible. The main text is therefore accompanied by many exercises and examples that 
have been worked out in great detail. This should make the book useful also for students 
wishing to study the subject on their own. 

Beginning the education in theoretical physics at the first university semester, and not as 
dictated by tradition after the first one and a half years in the third or fourth semester, has 
brought along quite a few changes as compared to the traditional courses in that discipline. 
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Especially necessary is a greater amalgamation between the actual physical problems and 
the necessary mathematics. Therefore, we treat in the first semester vector algebra and 
analysis, the solution of ordinary, linear differential equations, Newton’s mechanics of a 
mass point culminating in the discussion of Kepler’s laws (planetary motion), elements 
of astronomy, addressing modern research issues like the dark matter problem, and the 
mathematically simple mechanics of special relativity. 

Many explicitly worked-out examples and exercises illustrate the new concepts and 
methods and deepen the interrelationship between physics and mathematics. As a matter of 
fact, this first-semester course in theoretical mechanics is a precursor to theoretical physics. 
This changes significantly the content of the lectures of the second semester addressed in 
the volume Classical Mechanics: System of Particles and Hamiltonian Dynamics. 

The new mathematical tools are explained and exercised in many physical examples. In 
the lecturing praxis, the deepening of the exhibited material is carried out in a three-hour- 
per-week theoretica, that is, group exercises where eight or ten students solve the given 
exercises under the guidance of a tutor. 

Biographical and historical footnotes anchor the scientific development within the general 
context of scientific progress and evolution. In this context, I thank the publishers Harri 
Deutsch and F. A. Brockhaus (Brockhaus Enzyklopddie , F.A. Brockhaus, Wiesbaden — 
marked by [BR]) for giving permission to extract the biographical data of physicists and 
mathematicians from their publications. 

We should also mention that in preparing some early sections and exercises of our 
lectures we relied on the book Theory and Problems of Theoretical Mechanics, by Murray 
R. Spiegel, McGraw-Hill, New York, 1967. 

Over the years, we enjoyed the help of several students and collaborators, in particular, 
H. Angermuller, R Bergmann, H. Betz, W. Betz, G. Binnig (Nobel prize 1986), J. Briechle, 
M. Bundschuh, W. Caspar, C. v. Charewski, J. v. Czarnecki, R. Fickler, R. Fiedler, B. Fricke 
(now professor at Kassel University), C. Greiner (now professor at JWG-University, Frank- 
furt am Main), M. Greiner, W. Grosch, R. Heuer, E. Hoffmann, L. Kohaupt, N. Krug, 
P. Kurowski, H. Leber, H. J. Lustig, A. Mahn, B. Moreth, R. Morschel, B. Muller (now 
professor at Duke University, Durham, N.C.), H. Muller, H. Peitz, J. Rafelski (now pro- 
fessor at University of Arizona, Tuscon), G. Plunien, J. Reinhardt, M. Rufa, H. Schaller, 
D. Schebesta, H. J. Scheefer, H. Schwerin, M. Seiwert, G. Soff (now professor at Technical 
University Dresden), M. Soffel (now professor at Technical University Dresden), E. Stein 
(now professor at Maharishi University, Vlodrop, Netherlands), K. E. Stiebig, E. Stammler, 
H. Stock, H. Stormer (Nobel prize 1998), J. Wagner, and R. Zimmermann. They all made 
their way in science and society, and meanwhile work as professors at universities, as 
leaders in industry, and in other places. We particularly acknowledge the recent help of 
Dr. Sven Soff and Dr. Stefan Scherer during the preparation of the English manuscript. The 
figures were drawn by Mrs. A. Steidl. 

The English manuscript was copy-edited by Kristen Cassereau and the production of the 
book was supervised by Timothy Taylor of Springer- Verlag New York, Inc. 

Walter Greiner 

Johann Wolfgang Goethe-Universitat 

Frankfurt am Main, Germany 
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VECTOR CALCULUS 




1 Introduction and 
Basic Definitions 



Physical quantities that are completely determined by the specification of one numerical 
value and a unit are called 

scalars (e.g., mass, temperature, energy, wavelength). 

Quantities that for a complete description besides the numerical 
unit still need the specification of their direction are called 

vectors (e.g., force, velocity, acceleration, torque). 

A vector may be represented geometrically by an oriented 
distance, i.e., by a distance associated with a direction, such that 
holds; for example: Let A be the initial point and B the endpoint 
of the vector a (compare figure). 

The magnitude of the vector is then represented by the length of the distance AB. A 
vector is frequently described symbolically by a Latin letter with a small arrow attached to 
elucidate the vector character. Other possible representations make use of German letters 
or emphasize the quantity by bold printing. 

The magnitude of a vector a is written as: |a| = a. 

Definition: Two vectors a and b are called equal if 

1. Ia| = |b|, 

2. a ff b (aligned; parallel). 

We then write a = b. 

That means: All distances of equal length and equal The vectors a and b are equal, 
orientation are representations of the same vector on equal 
footing. Hence, the specific location of the vector in space is being disregarded. 

A vector with opposite direction but equal magnitude of a is denoted as —a. Oppositely 
equal vectors have the same length (|a| = | — a|) and are located on parallel straight lines 
but have opposite orientations; that is, they are antiparallel (a —a). If, for instance, 

a = AB, then —a = BA. 





Vector a pointing from A 
to B. 
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Addition: If two vectors a and b are added, the initial 

point of the one vector is brought by a parallel shift to 
coincide with the endpoint of the other one. The sum a+b. 
also called the resultant, then corresponds to the distance 
from the initial point of the first vector to the endpoint 
of the second one. This sum may also be found as the 
diagonal of the parallelogram formed by a and b (compare 
the figure). 

Rules of calculation: There hold 




Addition of the vectors a and b. 



a + b = b + a (commutation law) 

and 

(a + b) + c = a + (b + c) (association law), 
as is seen immediately (compare the figures). 




Illustration of the commutativity of the addi- Illustration of the associativity of the addi- 
tion of vectors. tion of vectors. 



Subtraction: The difference of two vectors a and b is defined as 

a — b = a + (— b). 

Zero (Null) vector: The vector difference a — a is denoted _ a 

as zero vector (or null vector): 

a 

a — a = 0ora — a = 0. 

The zero vector. 

The zero vector has magnitude 0; it is orientationless. 



Multiplication of a vector by a scalar: The product pa of a vector a by a scalar p, 

where p is a real number, is understood as the vector having the same orientation as a and 
the magnitude \pa\ = |p| • |a|. 
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Rules of calculation: 

q(pa) = p(qa) = qpa 
(p + q) a — pa + qa, 
p(a + b) — pa + pb. 



(where p and q are real). 





These rules are immediately intelligible 
and don’t need any further explanation. 



The multiplication of a vector a by a scalar p (in 
this case, p = 3). 



2 The Scalar Product 



The physical quantities force and path are oriented quantities and are represented by the 
vectors F and s. The mechanical work W performed by a force F along a straight path s is 

W — Fs cos (p = |F| | s | cos <p, 

where (p is the angle enclosed by F and s. W by itself, although originating from two 
vectors, is a scalar quantity. With a view on physical applications of this kind, we therefore 
define: 

The scalar product a • b of two vectors is understood as 
a • b = |a| ■ |b| ■ cos <p, 

where <p is the angle enclosed by a and b. a • b is a real number. Expressed by words, the 
scalar product is defined as follows: a • b = |a| multiplied by the projection of b onto a, or 
vice versa. 





Illustration of the scalar product. 



The visual meaning of the scalar product: 

magnitude of the projection of b onto a multiplied by |a|, or 
magnitude of the projection of a onto b multiplied by |b|. 
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THE SCALAR PRODUCT 2 



Properties of the scalar product: a • b takes its maximum value for q> equal to zero 

(cosO = 1, a parallel to b) 

a ■ b = |a| ■ |b|. 

For (p = 7t the scalar product takes its minimum value (cos7T = — 1, a antiparallel to b), 
namely 

a • b = — |a| • |b|. 

For <p — jr/2, a • b = 0 holds, even if a and b are nonzero (cos7r/2 = 0, 
a perpendicular to b); thus 

a • b = 0 if a_Lb. 

Rules of calculation: The following are true: 

a ■ b = b • a (commutativity); 

a (b + c) = a • b + a c (distributivity); 

p{ b ■ c) = (p b) ■ c (associativity). 

The first and last rules are immediately intelligible; the second rule is illustrated in the 
figure below. 

If b, c, a are not coplanar, the rule of distributivity may easily be visualized by a triangle 
located in space. The vector a may easily be visualized by a pencil or a pointing rod 
(compare the figures!). 



Unit vectors: Unit vectors are understood as vectors of magnitude 1. If a ^ 0, then 



is a unit vector pointing along a. Actually, the magnitude of e equals 1 since |e| = 
|a/|a|| = |a|/|a| = 1. A possibility frequently used in physics is to assign a direction to 




Illustration of the distributivity law. 



Illustration of the distributivity law in 
space. 
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a scalarly formulated equation by the unit vector. For example, the gravitational force has 
the magnitude 



mM 




It is acting along the connecting line between the two masses M 
and m, hence 




F is the force applied by the mass M to the mass m . Its direction 
is given by —e r = — r/|r|. Hence it is acting toward the mass M. 




The unit vector pointing 
from the big mass to the 
small mass is e r = r/|r|. 



Cartesian unit vectors: The unit vectors pointing along the positive x-, y-, and L-axes 

of a Cartesian coordinate frame are defined as follows: 

ei (in x -direction) or also i; 
e 2 (in y-direction) or also j; 
e 3 (in ’-direction) or also k. 

There exist two kinds of Cartesian coordinate frames, namely right-handed frames and 
left-handed frames (compare the figures below). 




right-handed system: k points into the di- 
rection of a right-handed screw when i i-x j 
is rotated along the shortest possible way. 



left-handed system: k points into the direc- 
tion of a left-handed screw when i i-v j is 
rotated along the shortest possible way. 




We shall always use only right-handed frames in these lectures! 

Orthonormality relations: i. j, k or ei , e 2 , e 3 will be used in the following always con- 

currently, depending on convenience. 

We now consider the properties of the Cartesian unit vectors with respect to formation of 
scalar products: Since the enclosed angle is each a right one, the following relations hold: 

i-i = j- j = k- k= 1 (because of <p = 0, hence cos 0=1); 
i-j = i- k = j- k= 0 (because of <p = n /2, hence cos7r/2 = 0). 
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These relations are combined by defining 



0 for v ^ fi, 

e /t • e„ = 8^ v , where 8 tlv = l 

[1 for v = ii, 

and is called the Kronecker symbol. 1 For the three-dimensional space, /i and v are running 
from 1 to 3, ei = i, ej = j, e 3 = k. 



1 Leopold Kronecker, b. Dec. 7, 1823. Liegnitz (Legnica) — d. Dec. 29, 1891, Berlin. Kronecker was a rich 
private person who moved to Berlin in 1855. He taught for many years at the university there, without having 
a chair. Only in 1883, after retirement of his teacher and friend Kummer, he took a professorship. His most 
important publications concern arithmetics, theory of ideals, number theory, and elliptic functions. Kronecker 
was the leading representative of the Berlin School, which claimed the necessity of arithmetization of the entire 
mathematics. 




3 Component 
Representation 
of a Vector 



The vector a, which is uniquely represented by the sum of 
vectors — in our example by the sum of the vectors b, c, d, f — 
is called the linear combination of the vectors (e.g., b, c, d, 
and f). The term “vectors” and their “linear combination” thus 
graphically form a closed polygon, the vector polygon. One 
may, of course, conclude from given vectors b, c, d on the 
linear combination that yields the arbitrary (but fixed) vector a. 

According to the definition introduced above, the vector a 
then must be a linear combination of the vectors b, c, d; 
thus 

a = q\\) + q 2 c + <? 3 d. 

qi,q 2 , and q 2 are denoted as components of the vector a with respect to b, c, d. The vectors 
b, c, d must be linearly independent, that is, none of the three vectors may be represented 
by the other two vectors. Otherwise not every arbitrary vector a could be combined out 
of the three basic vectors b, c, d. If, for example, d could be expressed by b and c, hence 
d = ab + /3c, then a = (<71 + q 2 a)b + (q 2 + q 2 f)c would always be confined to lie in the 
plane spanned by b and c. But an arbitrary vector a in general does not lie in this plane 
(e.g., points out of this plane). One says: The base b, c is incomplete for arbitrary vectors 
a. In the three-dimensional space one therefore always needs three basic vectors that are 
linearly independent (i.e., cannot be expressed by each other). 

Component representation of a vector in Cartesian coordinates: Any vector of the 

three-dimensional space may be represented as a linear combination of the Cartesian unit 




The vector polygon. 
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COMPONENT REPRESENTATION OF A VECTOR 3 



vectors i, j, k. This representation leads to 
simple and transparent calculations, due to 
the orthogonality relations. One then has 

a = a x i + ci y j + a-k, 

where a x = a ■ i, a y = a j. and a z = a-k 
are the projections of a onto the axes of the 
frame. The unit vectors i, j, k (or ej , e2, e^) 
are also called base vectors. 

Besides the representation as a sum of vec- 
tors along the unit vectors, the vector a still 
may be represented as 



a = (a x , a y , a z ) 

( a'\ 



a = 



a y 

\ a ~J 



(row notation), 



(column notation). 




The components of a vector are obtained by 
parallel projection. 



If the base vectors are known, it is sufficient to know the three components. 



Calculation of the magnitude of a vector from the components: According to the theo- 

rem of Pythagoras, the magnitude of a vector a is calculated from its Cartesian components 
as follows: 

|a| = + a 2 y + a\. 

Addition of vectors expressed by components: One has 

3 3 3 

a + h = y" a,ei + ^ Z?,e; = ^(a, + &,)e, 

1 = 1 1 = 1 1 = 1 

= («i + b\)t\ + (a2 + bi)^! + («3 + bi)^ 

= (fli + b \ , a2 + b2, «3 + b^) . 

Here both commutativity as well as associativity of vector addition have been used 
repeatedly. Thus, the components of the sum vector are the sums of the corresponding 
components of the individual vectors. 

The scalar product in component representation: One has 

a b = (a x \ + a y j + a-k) ■ (b x i + b y j + A-k). 

— a x b x i ■ i T a , b , l • j -|- a, b.\ • k - a ■ b , j • i -j- a y b y j • j T a y b z j • k 
+ a z b x k • i + a-b y k ■ j + a z b z k • k . 
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Taking into account the orthonormality relations (2. 1 ), we then get 

a • b = a x b x + a v b v + a-b z . (3.1) 

Finally, setting for the indices x = 1, for y = 2, and for z = 3, then one can write 

3 

a • b = y ^ ajbj. (3.2) 

i = l 

Hence, the scalar product of two vectors may be evaluated simply by multiplying the cor- 
responding components of the vectors by each other and summing over the three products. 



Problem 3.1 : Addition and subtraction of vectors 

A DC- 10 “flies” north-west at 930 km/h relative to ground. A strong breeze blows from the west with 
120 km/h relative to ground. 

What are the velocity and direction of flight of the plane, assuming that there is no wind deflection? 




The relative directions of wind and airplane velocity. 



Solution Let 

|v m | = 930 km/h, the velocity of the plane in the wind, 
| Vo | = the velocity of the plane without wind, 

| w| = the wind velocity. 

Now we can write 



w = 1 20 e x , 

v,„ = — 930cos(45°)e. v + 930sin(45°)e > , 
= — 657.61e x + 657. 61e,,, 

Vo = v m — w = — 777.61e v + 657.61e v 



v 0 = |v 0 | = 1018.39 km/h, 
I Tod 



tan^> = L—± — 0.846 
I Tod 



=>■ <p = 40.2°. 
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The position vector: A point P in space may be 

uniquely fixed by specifying the vector beginning at 
the origin of the coordinate frame and pointing to 
the point P as endpoint. 




O Z 



P(x,y,z ) 



The components of this vector, the position vec- 
tor, then correspond to the coordinates (x, y, z) of 
the point P . Thus, for the position vector, which is 
mostly abbreviated by r, there holds 



r = xi + yj + zk, or: r = (x, y, z); 



|r| = V * 2 + y 2 + z 2 . 



x 



The angle between two vectors : From the knowl- 

edge of the two possibilities for represent- The position vector and its coordinates. 

ing the scalar product 

a b = |a| |b| cos <p — a x b x + a y b y + a z b z , 
one obtains the following relation for the angle enclosed by a and b: 



a b 



a x b x + a y b y + a z b z 







4 The Vector Product 
(Axial Vector) 



One may define a further product between vectors. Here a new vector arises that is defined 
as follows. 



Definition: The vector product of two vectors a and b is the vector 

a x b = (|a| • |b| sin<p)n. 



( 4 . 1 ) 



where n is the unit vector being perpendicular to the plane fixed by a and b, and pointing 
out of the plane as a right-handed helix when rotating the first vector of the product into 
the second vector. Note that the rotation has to be per- 
formed along the shortest path. 

The magnitude of the vector product is equal to the area 
of the parallelogram spanned by a and b, as is seen from 
the figure. 




F = | a x b| = | a | ■ |b| sin^> = ab sirup, 



Geometrical interpretation of the 
absolute value of the vector prod- 
uct as area. 



Properties of the vector product: a x b takes its maximum magnitude for <p = jr/2, 

sin(7r/2) = 1, a perpendicular to b, |a x b| = |a| |b|. 
a x b vanishes for ip — 0 (sin 0 = 0, a parallel to b). 

I if a b or means antiparallel) 

if a ft b. (ff means parallel) 

The formula also includes the special case a = b, thus 

a x a = 0. 
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THE VECTOR PRODUCT (AXIAL VECTOR) 4 



Notations: 

O represents a vector perpendicular to the drawing plane and pointing out of the plane 
(arrowhead). 

® represents a vector perpendicular to the drawing plane and pointing into the plane 
(arrowbase). 

Rules of calculation: The vector product has the following properties: 

I. axb = — b x a (no commutativity); 

II. a x (b + c) = a x b + a x c (distributivity); ^ 

III. a x (b x c) (a x b) x c (no associativity); 

IV. p( a xb) = (pa) xb = ax (pb). 

Rule I follows immediately from the definition of the vector product (compare with the 
figure). 

Rule III: The vector on the left side lies in the plane spanned by the vectors b and c; the 
vector on the right side is in the plane spanned by a and b. The subsequent example also 
shows that associativity does not hold. One has ei x (e 2 x ej) — 0, but (ei x ea) x ej = — ei . 
Rule II: The proof is given in two steps: 

1 . Let a be perpendicular (_L) on b and c, that is, a • b = a ■ c = 0. Then a x (b + c) = 
a x b + a x c. The proof for that may be read off immediately from the two figures, 
a x b stands _L on b and a, is rotated against b by 90°, and is longer than b by the factor 
|a|. The situation is similar for axe and ax (b + c). The parallelogram of the vectors 
a x b, a x c, a x (b + c) emerges from that of the vectors b, c, (b + c) by a rotation 
about a by 90° and subsequent stretching by |a|. 




Illustration of the calculational rule I. The direction of rotation is shown. 
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Plan view from top of the special case: alb 
and ale. 



2. We now decompose in the general case 

b = bj_ + b||, 

C = C ± + C||, 

that is, b and c into components _L and 1 1 to a (compare 
with the figure). 

Then, on the one hand, the following holds: 




axb = (|a| • |b| ■ simp) 



a x b 

|a x b| ’ 



and, on the other hand, 



The general case: the vectors 
b and c are both decomposed 
into components parallel (|| ) and 
perpendicular (_L) to a. 



a x b 

a x b± = (|a| • |b x |)- — 

|a x b| 

a x b 

= (|a| • |b| • sin(c) - 

|a x b 

and therefore 



a x b = a x b ± . (4.3) 

This holds for any arbitrary vector b. Therefore, one immediately concludes a x c = 
a x c_l . We may then conclude 

a x (b + c) = a x (b + c)x = a x (bj_ + c±) 

= a x b± + axci (because of the special case in 1) 

= axb + axc (because of (4.3)) 

q.e.d. 
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Rule IV: The rule for multiplication by a scalar p is immediately evident if we remind 
ourselves of the meaning of pa. 

Vector products of the Cartesian unit vectors: There holds 

ixi = jxj = kxk = 0, and i x j = k, 
or ei x ei = e2 x e2 = e3 x e3 = 0, and ei x ei = e3. (4.4) 

This product satisfies the cyclic permutability. For an anticyclic permutation one has to 
multiply by the factor — 1, for example, j x i = — k. 

Vector product in components: We now denote the Cartesian unit vectors by ei, e2, e3 

instead of i, j, k. 

Let 

3 3 

a = aid + + #3^3 = a,e, and b = b,t\ 

i=i i=i 

be two arbitrary vectors. When forming the vector product of the two vectors a = Y^=i a > e > 
and b = Y^=i one obtains 

a x b = (aiei + 02^2 + a 3 e 3 ) x (b iei + b 3 ^2 + b 363) 

= < 21^263 — «2^i e 3 + «2^3 e i — a 3 b 3 ^i + aj,b\^2 ~ a\b 3 £2 (4-5) 

= (a 2 b 3 - a 3 b 2 )ei + (a 3 bi - a\b 3 )e 2 + (a\b 2 - a2bi)e 3 . 

It is now practical to introduce the determinant notation. 



Determinants: A rectangular array of numbers is called a matrix (see the figure). 



/ «11 


column 

i 

fll2 ■ ■ • 


■ a \q \ 




«21 


fl22 • • • 


• &2q 


row 


V Gp\ 


fl/j2 • ■ • 


a pq / 





For the case q — p, the matrix is called quadratic. One then can assign a numerical value 
D to it, called a determinant. It is defined as follows: 



I. 

II. 



det(fln) = |an| = an; 



/an 


fll2\ _ 


ail 


fll2 


\fl2l 


022/ 


fl21 


«22 



— a 1 1 ^22 — a i2«2L 



(4.6) 
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an 


Cl\2 


an 




a n 


an 


013 


an 


an 


an 


= 


021 


«22 


023 


y«31 


an 


an ) 




«31 


«32 


033 





022 


023 




021 


023 


+ O 13 


021 


022 


flu 


O 32 


033 


— fln 


031 


033 


«31 


O 32 





The evaluation of the 3x3 determinants may be simplified by using the so-called Sarrus 
rule. 1 The procedure is: Establish an additional auxiliary matrix by writing the first two 
columns of the original matrix once again to its right side, and form the product terms, 
involving signs, according to the following scheme. 




Multiple-row determinants may be reduced to determinants of lower order by expansion 
with respect to a row or column (formation of subdeterminants), analogous to (eq. (4.6), 
III). We will see this method at work in Example 4.4 on the Laplace expansion theorem. 

Rules of calculation: The most important rules for calculations involving determi- 

nants are 

1. If two rows or columns of the quadratic matrix are identical or proportional to each 
other, then the determinant of this matrix = 0. 

2. When permuting any two neighboring rows or columns, the sign of the determinant 
changes. 

3. The determinant of the matrix reflected at the main diagonal (also called the transposed 
matrix) is equal to the original determinant. 



1 Pierre Frederic Sarrus, b. 1798 — d. 1861, Saint Affriques. Sarrus was professor of mathematics in Strasbourg 
from 1826 untill 1856. He dealt mainly with the numerical solution of equations with several unknowns (1832), 
with multiple integrals ( 1 842), and with the determination of comet orbits (1843). The rule for evaluating three-row 
determinants is named after him. 
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4. The expansion theorem that has been used in equation (4.6), III, with respect to the first 
row, holds for the first column in the same way. 

These rules may easily be checked explicitly in the cases quoted above (eq. (4.6), I, II, 
III). The cases i -m are the most important ones in the present context. The properties 
of the 3x3 determinants will be outlined and discussed in more detail in the context of 
Problem 4.3. The rules hold, however, in general for arbitrary determinants. 

The vector product (4.5) may now be written as a three-row determinant: 



ei e 2 e 3 

a x b = |«i a 2 «3 
bi b 2 b 3 

= efa 2 b 3 - a 3 b 2 ) + e 2 (a 3 bi - a\b 3 ) + e 3 (a 1 b 2 - a 2 b 0 . 



( 4 . 7 ) 



If the two vectors of the cross product are equal, then the two lower rows of the determinant 
are also equal, and the vector product vanishes. 

Further, one may easily check based on equation (4.6), III, that the sign of the determinant 
changes under permutation of rows (or columns). This corresponds to the anticommutativity 
of the vector product. 



Representation of the product vector: As we already stated in the definition of the 

vector product, the magnitude of the product vector may be visualized by a distance but 
better by the area of the parallelogram formed by the vectors. This vector is not determined 
by its length and orientation only (such vectors are called polar vectors) but is called 
an axial vector. To understand this difference, we consider space reflections: We thereby 
change from the components a i, a 2 , a 3 to the new base vectors a\ — —a 3 , al, = —a 2 , 
flj = —c/ 3 . The vector a is thus reflected at the origin. Under a space reflection, which 
is also called a parity transformation , a polar vector changes its sign: a — > —a. An axial 
vector, on the contrary, remains unchanged: axb = (—a) x (— b). 

The invention of these new vectors is necessary since certain physical quantities need a 
handedness for a complete description. The handedness is taken into account by an axial 
vector. Such kinds of quantities are, for instance, the angular velocity and the angular 
momentum. One should get straight in one’s mind that a handedness remains unchanged 
under a space reflection! 

An axial vector may, however, also be represented by an oriented distance. 

The double-vector product: The vector product a x (b x c) is called the double-vector 

product. To evaluate it, we denote the components of b x c as follows: Let 

(b x c) A be the .v -component, 

(b x c) y be the y -component, and 
(b x c) z be the ^-component. 
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For the .r-component of the double-vector product, it then follows that 

(a x (b x c)). v = a y (b x c) z -a,(bx c) y 

Cly{b x Cy byC X ) d- (ll-Cy l) X C- ) . 

We add a x b x c x — a x b x c x = 0 and obtain 

(a x (b x c)). v = b x (a x c x + a y c y + a z c z ) - c x (a x b x + a y b y + a z b z ) 

— b x (a ■ c) — c v (a • b). 

Analogous considerations for the y- and ^-components of a x (b x c) yield the 

Graliniann expansion theorem: One has 

a x (b x c) = (a ■ c)b — (a • b)c, 
while 

(a x b) x c = (a ■ c)b — (b • c)a. (4.8) 

This is another proof of the fact that the vector product is not associative (see (4.2), III). 



Problem 4.1 : Vector product 

(a) The vector (1, a, b ) is perpendicular to the two vectors (4, 3, 0) and (5, 1, 7). Find a and b. 

(b) Evaluate in Cartesian coordinates the vector product a x b for a = (1, 7, 0) and b = (1, 1, 1). 

(c) Show that 

(a x b) 2 = a 2 b 2 - (a • b) 2 . 

Solution (a) It must hold that (1 , a, b) ■ (4, 3,0) = 0 and (1, a,b) ■ ( 5, 1, 7) = 0. This yields the two equations 

4 11 

4 + 3a = 0 and 5 + a + 1b = 0 =>■ a = ,£>= . 

3 21 

(b) 

(a x b), = (a y b z - a z b y ) = 7: 

(a x b) y = (a z b x - a x b z ) = -1; 

(a x b) z = (a x b y - a y b x ) = -6. 

(c) 

(a x b) 2 = (|a| • |b| • sin<p ■ e„) 2 = |a| 2 |b| 2 sin 2 (p(e„) 2 
= |a| 2 |b| 2 (l - cos 2 <p)(e„) 2 
= a 2 b 2 — (a • b) 2 



Here (p (a, b) and e„ is the unit vector along a x b. 
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Problem 4 . 2 : Proof of theorems on determinants 

The most important theorems on determinants are as follows: 

(a) Under permutation of rows and columns (reflection at the main diagonal), the value of a deter- 
minant remains unchanged. 

(b) Under permutation of two arbitrary neighboring rows, the sign of the determinant changes. 

(c) If all elements of a row contain a common factor c, then it may be pulled out of the determinant. 

(d) If two rows of a determinant are proportional to each other, then the determinant = 0. 

(e) The value of a determinant remains unchanged when adding a multiple of any row to another 
row. 

Check these rules for a general 3x3 determinant. 

Solution From the lecture we know the definition of the 3-determinant: 





Oil 


012 


013 


D = 


021 


022 


023 




031 


032 


033 



= All (<322^33 — 023032) ~ 012)021033 — 023^31) + 013(021032 — O22O31) 

— O11O22O33 — O11O23O32 — O12O21O33 + 012023031 + O13O21O32 — 013022031. ( 4 . 9 ) 

(a) Permutation of rows and columns of D (reflection at the main diagonal) leads to 





Oil 


O21 


O31 


D = 


O12 


O22 


O32 




O13 


O23 


O33 



= 011(022033 — 032023) — 021(012033 — 032013) + 031(012023 — O22O13) 

— 011022033 — 0iifl32023 — 031012033 + 021033013 + 031012033 — O31O22O13 

— 011(022033 — 023032) — 012(021033 — O23O31) + 013(021032 — O22O31). 

A comparison with D (see above) yields 

D=D. ( 4 . 10 ) 

(b) Permutation of, for example, the second and third rows of D yields 





Oil 


012 


013 


D' = 


031 


032 


033 




021 


022 


023 



= On (032023 — O33O22) — 012(031023 — 033021) + 013(031022 — 033021) 

— O11O32O23 — Onfl33022 — O12O31O23 + O12O33O21 + 0l3O3iO22 — 013032021- 
By means of 4.9, one immediately concludes that 



D' = -D. 



( 4 . 11 ) 
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This means: When we permute the second and third rows, the determinant changes its sign. Similarly, 
one may check that for permuting other rows. From (a) the same result follows for the columns: 
When we permute neighboring columns, the determinant also changes its sign. 

(c) We investigate 



Oil 


012 


013 


C021 


C022 


C023 


031 


a 32 


033 



= an (0222033 — Cfl23^32) — 0 l 2 (CO 2 lO 33 — CO 23 O 31 ) + Ol3(c02l032 — C 02203 l) 

= c [ail(a 22033 — 023032 ) — 012(021033 — O 23 O 31 ) + 013(021032 — O 22 O 31 )] 

and compare with 4.9. Obviously, 

D" = cD. (4.12) 

(d) For example, let the third row be proportional to the second row; thus 



Oil 


a 12 


013 


021 


022 


023 


/.02 1 


A022 


kfl 23 



= On (A.U22O23 — /-O03O22) — o 1 2 (/-O21 023 — /-(223O21 ) + o 11 (/.i/ 2 i O22 — /.O22O2] ) 

= 0. (4.13) 

Similarly, one may check the assertion for the proportionality of other rows. From (a), it follows 
immediately that the determinant also vanishes if two columns are proportional to each other. 

(e) We add. for example, a multiple of the first row to the second row. Then 





Oil O12 O13 




5 " = 


O21 4 “ A.OH O22 4 ~ i.o ] 2 O23 4 ~ /.3i ] 3 






031 032 033 




— fill ^(022 + ^012)033 — (023 + A.fli 3 )a 32 


] 




— 012 1 


^(021 4 - kfln)a 33 — (023 4 - kai 3 )fl 3 i j 




4 - 013 1 


^(021 4 - Xfln)a 32 — (022 4 - ^012)031 j 



— O11O22O33 + ^'• a ll fl 12tf33 — 011023032 ~ ^ a l 1O13O32 

— O12O21O33 — ^-^ 12 ^ 11^33 4 “ O12O23O31 4 “ ^-0 I2O13O31 

+ fll3021<t32 + ^ai3aua32 ~ 013022031 — A.£Zi3«i2<5t31 
= «1 1 (C(22«33 — fl23<t32) — 012(021033 — 023031) + 013(021032 — fl2203l). 

A comparison with 4.9 yields the assertion 



D" = D. 



(4.14) 
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Problem 4.3: Determinants 

Calculate using the theorems on determinants: 





X 


x + 1 


x + 2 


a 


d 


xa + yd 




4 


5 


22 


(a) 


0 


1 


2 (b) 


b 


e 


xb + ye 


(c) 


8 


11 


44 




3 


3 


3 


c 


f 


xc + yf 




3 


7 


1 



Solution (a) We form the linear combination 

x 

a ■ (2. row) + ■ (3. row) with a = 1, /S = — 

and obtain (,r, x + 1, x + 2), thus just the first row. From (i) and (ii) of (4.6) it follows that the 
determinant is always equal to zero. 

(b) The third column is a linear combination of the first and second columns with the factors x and y : 

( a\ / d \ / xa + yd \ 

From this it follows that the determinant becomes zero. 

(c) We expand with respect to the first row: 

11 44 8 44 8 11 

4 - 5 +22 = 4( — 297) — 5( — 124) + 22(23) = —62. 

7 1 3 1 3 7 

Example 4.4: Laplace expansion theorem 

Let A — ( da ) bean x n matrix, and So, be the submatrices of A obtained by erasing the ; th row and 
the kth column of the matrix A. The matrices S,a thus are (n — 1) x (n — 1) matrices. For each i with 
1 < ( < n, it holds that 

n 

det A = T>l) l+ % det 5,* (expansion with respect to (th row) 

k= t 

and also 

n 

det A = ^^(— l) l+k a k[ det 5 a,- (expansion with respect to i th column). 

A=1 

We check the theorem explicitly for 3-determinants and expand at first the general 3x3 determinant: 
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an «i2 flu 

Expansion of det A — a 2 1 022 <223 with respect to the first row yields 

(231 <232 «33 

det A = (— l) 1 + 1 aii 5 n + (— l) 1 + "ai2‘Si2 + (— D 1 + 3 «i 3 ■S'i 3 





O 22 


a 23 




021 


023 


+ «13 


02 1 


a 22 


On 






— a l2 












a 32 


(?33 




Ojl 


a 33 




O 3 I 


a 32 



(4.15) 

(4.16) 



Expansion of the 3-determinant with respect to the second column yields 

det A = (— l) 1+ “(2l2Sl2 + ( — 1 ) -+- «22*S , 22 + (— l) 3+ “fl32>S , 32. (4-17) 

The first term on the right side is identical with the second term of 4.15. The last two terms of 4.17 
read explicitly 




Obviously, 4.18 and 4.19 coincide. Hence, it is clear that the expansions of the 3-determinant with 
respect to the first row and the second column, respectively, yield the same. Similarly, one may 
verify that the expansion with respect to other rows or columns leads to the same result. Hence, the 
expansion theorem for 3-determinants is seen to be valid. 

We still evaluate the 3x3 determinant by expanding with respect to the second row, and subse- 
quently with respect to the second column, for the example of the determinant 

4 5 22 

det A = 8 11 44 

3 7 1 

This yields 

(a) Expansion with respect to the second row: 

det A = (— 1) 2+I a2l5 , 21 + (— 1) -+- <222*S22 + (— l) -+3 fl23^23 
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(4.20) 
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(b) Expansion with respect to the second column: 

det A = (— l) 2+1 «l 2 Sl 2 + ( — 1 ) -+ “<322‘S , 22 + (— 1)- +3 032 S32 





021 


023 


+ «22 


Oil 


Ol3 




Oil 


013 


o 12 










— O 32 








«31 


033 




031 


033 




021 


023 





8 44 




4 22 
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OO 



( 4 . 21 ) 




5 The Triple Scalar 
Product 



Definition: The triple scalar product of the three vectors a, b, and c is defined as 

a • (b x c) , 

that is, a combination of a scalar and vector product. The triple scalar product is therefore 
also denoted as a mixed product. The triple scalar product is a scalar. 



Triple scalar product in component representation: 

a • (b x c) = (ai, a 2 , a 3 ) ■ [(£>1, b 2 , b 3 ) x (ci, c 2 , c 3 )] 



= (ai,a 2 , a 3 ) 



ei e 2 e 3 

b\ b 2 b 3 

Ci c 2 c 3 



= (at, a 2 , a 3 ) ■ ( b 2 c 3 - b 3 c 2 , -b\c 3 + b 3 c i, £>ic 2 - b 2 c{) 
= a.\(b 2 c 3 - b 3 c 2 ) - a 2 (b\C 3 - b 3 Ci) + a 3 (b x c 2 - b 2 Ci). 



The three terms may again be combined to a determinant, such that 



a • (b x c) = 



G\ Cl2 

b\ b 2 b 3 

ci c 2 c 3 



= (a x b) • c. 



(5.1) 



Cyclic permutability: The factors of the triple scalar product may be permuted cyclically. 

One has 

a • (b x c) = b • (c x a) = c • (a x b). 
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These rules may be confirmed easily by successive permutations of the rows in the deter- 
minant (5.1). The following simplified notation for the triple scalar product may be found 
occasionally in the literature: 

a ■ (b x c) = [abc] = [bca] = [cab]. 

Geometrically, the triple scalar product represents 
the volume 

V = a ■ (b x c) = a cos cp be sin y 
= abc cos (p siny 

of a parallelepipedon formed by the three vectors (see 
figure). 

Note: The volume has a positive sign (+) if a 

lies on the side of b x c, but a negative sign (— ) if a lies on the side of — b x c . Hence the 
volume might be associated with a sign. In general, however, this choice is not used, and a 
positive sign is always required. This is achieved by the definition V = a • (b x c)|. 




Properties of the triple scalar product: From 

TV 



a • (b x c) = 0 follows 




and / or y = 0, 



(5.2) 



that is, the three vectors are coplanar or (and) two vectors lie on a straight line. 

This is again a very clear proof of the theorems on determinants already mentioned 
above: 

1 . If two row vectors (or column vectors) are equal or proportional to each other, then the 
determinant equals zero. 

2. When we permute two neighboring rows, the determinant changes by a factor (—1). 



6 Application of Vector 
Calculus 



Application in mathematics: 



Problem 6.1 : Distance vector 

Calculate the length of the vector a that represents the distance vector 
between the points rj and r 2 . 

Solution a = r 2 - r t 

— fe e i + + £2^3) — (x 1 e ! + yie 2 + z^) 

= O2 - *t)ei + (y 2 - Vi)e 2 + (z 2 - zi)e 3 ; 
hence a reads in row notation 
a = (x 2 - xu y 2 - yi, z 2 - zi), 
and the magnitude of a is therefore 

|a| = -/ (x 2 - X \) 2 + (y 2 - yf ) 2 + (z 2 - Z1) 2 . 




The distance vector be- 
tween the points ri and 

t*2- 



Problem 6.2: Projection of a vector onto another vector 

Given 

a = (2, 1, 1), 
b = (1, -2, 2), 
c = (3, -4, 2), 

what is the absolute value of the projection of the sum (a + b) onto 
the vector c ? 




The projection of the sum 
a + b onto the vector c. 



Solution This projection is given by the scalar product of (a + b) and the unit vector e r along c. 

c (3, -4, 2) 

Cc _ Tci _ 
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(a + b) = (2+ 1, 1 -2, 1 + 2), 



(a + b) • e c = 



3 ■ 3 + (— 1) ■ (—4) + 3-2 
V29 



19 

V29' 



Problem 6.3: Equations of a straight line and of a plane 

Let the points A and B be given by their position vectors a and 
b. What is the equation of the straight line through A and B? 

Solution The straight line AB is parallel to (b — a). Moreover, it passes 
through point A. Hence, the equation determining any position 
vector x of a point X on the desired straight line reads 

x = a + r(b — a), 

with t being a real number (running parameter — oo < t < oo). 
If two points A and B are not given but one point A and a vector 
u specifying the orientation of the straight line are given, the 
equation of the straight line reads 




line. 



x = a + fu . 



This is called the point-direction form of the equation of a straight line. 
Example: 



a = (fli, t/2, 03), u = («i, u 2 , uf), 

X = (fli + til ] , @2 + lW2 , + tuf) 

= (x, y,z). 

A plane in space may be fixed by specifying be- 
sides the position vector a and the orientation vector 
u still a second orientation vector v: 

x E = a + mi + k\, 

where u v and also u f) v and k,t e R. The 
notation ff and indicates that u and v are neither 
parallel nor antiparallel. 

This is the point-direction form of the equation of 
the plane. 




Representation of a plane in space spanned 
by the vectors u and v from point P 0 . 



Example 6.4: The cosine theorem 

The cosine law of plane trigonometry is obtained by scalar multiplication 
of the equation c = a — b by itself: 

c • c = (a — b) ■ (a — b) = a 2 + b 2 — 2a ■ b 
= a 2 + b 2 — 2 ab cos y. 

=> c 2 = a 2 + b 2 — 2 ab cos y. 

For y = rr/2 there results the theorem of Pythagoras. 




a 



The vectors a, b, and c 
characterize the sides of 
the triangle. 
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Example 6.5: The theorem of Thales 

In order to prove the theorem of Thales 1 we introduce the 
following vectors according to the sketch: 

~M% = —MB = a, Wfc = b. 

It holds that 

|a| = |b|, BC = a + b, and a£ = b — a. 

The scalar product (a + b) • (a — b) has the value The theorem of Thales, demonstrated 

with the help of vectors. 

(a — b) • (a + b) = a 2 — b 2 = |a| 2 — |b| 2 = 0. 

For the angle enclosed by (a + b) and (a — b), it follows that i? = it / 2 or 
(a + b) ± (a — b) (theorem of Thales). 




Example 6.6: The rotation matrix 

The opposite figure shows into which vectors e'[ and e 2 the 
Cartesian unit vectors ei and e 2 are transformed under a 
rotation in the x, y-plane by the angle fi around the z-axis: 

e'j = ei cos fi + e 2 sin fi + e 3 • 0 

e 2 = ei (— sin fi) + e 2 cos /? + e 3 ■ 0 (6.1 ) 

= ei • 0 -P c 2 ■ 0 -P e 3 ■ 1 . 

This system of equations may be written in matrix form (see 
equation 6.7): 

( e'j \ /cos fi sin fi 0 \ / ei \ 

e 2 J = I — sin fi cos fi 0 ) • ( c 2 ) 




Case 1 : vector r stays at rest; the co- 
ordinate system is rotated. 



( r/nei + di 2 e 2 + ^1363 \ 

tf 2 1 ei + d 22 e 2 + rf 23 e 3 (6.2) 

^3 1 e i + ^32^2 + <^33 e 3 ' 



or briefly as 

3 

e v = 

M=1 



Earned after Thales ofMilet, b. about 624 BC — d. 546 BC. He is the first representative of the Ionic School. 
According to writings he did far travels (e.g., to Egypt) and was very active as a politician. The theorem named 
after him was for the first time strictly formulated by him. 



30 



APPLICATION OF VECTOR CALCULUS 6 



where 

( cos (8 
-sin/3 
0 

represent the direction cosines. Note that sirup = cos (cp — n /2). The matrix ( d vil ) describes the 
transformation of the base vectors. For a rotation in the three-dimensional space in the x , y-plane 
(i.e., about the z-axis), the rotation matrix reads 

( cos /3 sin yS 0 \ 

— sin/i cos fi 0 | ={d vil ). (6.3) 

0 0 1 / 



sin /3 0 

cos /3 0 

0 1 



Case 1: r = r' fixed in space. If r = r' is fixed in space but the coordinate frame rotates, one has 

V II 

Multiplication of this equation with e' ( isolates x ' : 

X fl = x v (e v • e (j ) = ^ d /iV x v . 

V V 

Thus, the transformation of the components of a position vector that is kept fixed in space is given by 
r' = Dr, (6.4) 

where D denotes the rotation matrix. Explicitly, this means because of x[ = x', x' n = y', x' 3 = z': 



x' 

y' 

z! 



( cos sin /i 0 

— sin yS cos y8 0 

0 0 1 




(6.5) 



The addendum “new base” at the column tuple shall indicate that the components x ' , y' , z! of the 
column tuple are to be interpreted as coefficients of the base vectors e' l5 e' 2 , and e' 3 . Written explicitly, 
the vector in the new basis thus reads 



r' = x'e'j + y'e 2 + z'e', . 

Case 2: r is tightly fixed to the rotating coordinate frame. 

Thus, r rotates with the coordinate frame. This means 

^^jr v e v = y^T a e,i 

V !L 

^ x ii = y x v (e^ ' Gyj) 

V 

= y d v/l x v 

V 

= V. d^yXy . ( 6 . 6 ) 




Case 2: vector r is rotated together 
with the coordinate system. 
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x'^ are the new components of the rotated vector with 
respect to the fixed system e M : x v are old components of 
the vector with respect to the fixed system e (1 . 

Note: Both x' v as well as x v in this case are defined in the old system (base e^). They denote the 

components of the new (rotated) and old (not rotated) vector, respectively! 

In the preceding we have already used the matrix multiplication. It shall once again be clearly 
defined here. 



Definition of the matrix product: The common element C, ( of the row i and the column j of the 

product matrix C = A ■ B is obtained by forming the sum 

C t] (6.7) 

k 

where A and B are the factor matrices. 

Thus, the components of a vector a = (ai, a 2 , 03 ) under rotations of the coordinate frame would 
change to 



= a = 



cos /I 


sin p 


°\ 


/a, 


— sin /S 


cos P 


0 • 


1 a 2 


0 


0 


l/ 


\Cl3 


cos pa [ + sin ft a 2 \ 




— sin p ai + cos 


Pa 2 I 


, 



Cl 3 



a t i — ^ 



Cly . 



The vector itself remains fixed in space. Its components change, however, because the base was 
rotated (case 1). If the vector would rotate (case 2), then we would obtain according to 6.6 



( ai \ 


1 


1 a ' 2 


= 


\ a ' 3 / 


\ 

new base 



cos f) a 1 — sin P a 2 ' 



a 3 






where d ^ = c/ V(1 is the transposed rotation matrix. The transposed of a matrix is simply the matrix 
reflected at the main diagonal (from the upper left to the lower right comer). 



Application in physics: 



Problem 6.7: Superposition of forces 

Four coplanar forces are acting at the point 0, as shown 
in the sketch. 

Calculate the net force acting at the point 0! 




Solution 



a = (-95.3, 55.0) N, b = (-150.4, -54.7) N, 
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c = (70.7, 70.7) N, d= (80.0, 0.0) N (n = Newton = 1 . 

It holds that 



F ges = a + b + c + d = (-95.0, 71.0) N, 
I F ges I = V95.0 2 + 71.0 2 N = 118. 6 N. 

We remember that 



a + b = ^2 a,e, + ^ (ye,- = ^(a, + b,)e, 



(fli +bi,a 2 + b 2 , a 3 + b 3 ). 



Graphical determination of the force: Representation by means 
of polygon of forces. 

The angle fi enclosed by F and the x-axis may be calculated 
easily. One has 

F v 71.0 

F = (—95.0, 71.0) N, — = tan /l = ; 

F x 95.0 

from there it follows that /3 = 143°. 




Graphical determination of the 
net force F. 



Example 6.8: Equilibrium condition for a rigid body without fixed rotational axis 

A rigid body is under the action of the forces F, at the positions r,- . We investigate the equilibrium at 
the point A (position vector a) the body may rotate about. All forces F,- are now added and subtracted 
at A such that nothing is changed in total. 





A rigid body is in equilibrium with respect to A pair of forces results in a torque M = 

point A if the sum of all torques with respect to r x F. These forces set the body on 

A and the sum of all forces in A vanish. If this which they act into rotation, 
condition is valid in A, it is also valid in every 
point B. 
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The procedure is illustrated for the force Fj in the figure on the next page; for the other forces 
we proceed in the same way. Now the forces Fi at ri and — Fi at a are forming a pair of forces that 
generates the torque (compare with the Problem 6.9) 



and will rotate the body. Similarly, all other forces F,- (at r, ) and — F,- (at a) are forming pairs of forces 
with the torques 



The question arises of whether an equilibrium at point A Equilibrium with respect to point A 
also means an equilibrium at point B. To answer it, we ' m P'' es equilibrium with respect to 
recalculate the eqs. 6.14, 6.15 to the point B (eqs. 6.12 P°' nt 
and 6.13). We realize: 6.14 is identical with 6.12. Further, it holds that 

M(b) = £>,- - b) x F, = £(r - (a + c)) x F, 

i i 

= y^(r, -a)xF;-^cx F, 

i i 

= M(a) - c x ^F,- = M(a) = 0. 



Mi (a) = (ri - a) x Fi 



(6.8) 



M,- (a) = (r,- - a) x F;. 

The total force acting at the point A is therefore 



(6.9) 




( 6 . 10 ) 



and the total torque about A is 




( 6 . 11 ) 



At the point B (position vector b ) a similar construction would yield the total force 




( 6 . 12 ) 



The total force F tries to accelerate the body as a whole. 
The total torque tries to rotate the body. If there shall be 
equilibrium with respect to point A (position vector a), then 
both the total force F and the total torque M(a) must vanish: 



and the total torque about B, 



F = 0, 
M(a) = 0. 




(6.13) 



(6.14) 

(6.15) 0 




=0 because of 6.14 



Therefore we may claim: If the equilibrium conditions 6.14,6.15 are fulfilled at a point A, then they 
also hold at any other point B. 
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Problem 6.9: Force and torque 



The following external forces are acting on a body: 

Fi = (10, 2, — 1) N at point P\ (2, 0, 0) cm, 
kg m 

(N = Newton = l-=— ) 

s 2 

F 2 = (0, 0, 5)N at point P 2 (l, 3, 0) cm 
and 

F 3 = (—6, 1, 8) N at point P 3 ( 6, 8, 1) cm. 



Find 




Illustration of the torque induced 
by two forces. 



(a) components, magnitude, and orientation of the resulting force F, 

(b) the torque with respect to P 2 . 



Remark (kp = kilopond). The kilopond is no longer a legal unit, all forces are measured in 
Newtons (N): 



m ,g cm , 

1 N = 1 kg • — = 10 5 — = 10 5 dyn, 

s 2 s 2 



1 kp = 9.81 N. 



Solution (a) 



F = Fj +F 2 + F 3 = (4, 3, 12) N, 

|F| = y/ 4 2 + 3 2 + 12 2 N = 13 N, 

cos Pi = — = 0.308, A =72°, 

|F| y 

Fy 

COS fa = — = 0.231, 62 = 11°, 

H |F| y 

cos p 3 = — = 0.923, fly = 23°. 
y |F| 

These are the direction cosines of the force. They describe the direction of force 



F 

n = — = (cos f)i, cos /3 2 , cos /3 3 ) = (0, 308; 0, 231; 0, 923). 



(b) The torque of a force F p acting at point 
P(x, y, z), that is at the position r = (x, y, z), is 
defined with respect to the coordinate origin (center 
of rotation) as the vector 

M = rxF p . 

Here r is the position vector from the cen- 
ter of rotation to the action point of the force 
F p . The magnitude of M is obviously given by 




Center of 
rotation 

The absolute value of the torque results from 
the force component perpendicular to the dis- 
tance vector. 
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M = r- F', where F' = F sin & is the force component perpendicular to the position vector (compare 
in the figure). This may also be expressed as follows: M = distance from center of rotation to action 
point of the force times force component perpendicular to the distance vector. 

This torque M is also caused by a pair of forces, as 
discussed in Example 6.8. If one adds at the center of 
rotation the forces — F and F, in total 0 (compare to the 
figure), then the forces — F at the center of rotation and F 
at r form a pair of forces. The force F acting at the center 
of rotation presses onto the bearing of the rotation axis and 
is received there. 

If several forces F v are acting on the rigid body at 
the points r„, the total torque is 

M = M„ = r„ x F v . 

V V 

In our example 

M = (it x Fi) + (r 2 x F 2 ) + (r 3 x F 3 ), 
r i = Pi - P 2 = (1, -3, 0) cm, 
r 2 = p 2 — p 2 = (0, 0, 0) cm, 
r 3 = p 3 - p 2 = (5, 5, 1) cm, 

where p; , p 2 , and p 3 are the position vectors of the points Pj , P 2 , and P 3 . Hence one obtains 

M| = it x Fj = (3, 1, 32) N cm, 

M 2 = r 2 x F 2 = ON cm, 

M 3 = r 3 x F 3 = (39, -46, 35) N cm. 

The total torque is 

M = M! + M 2 + M 3 = (42, -45, 67) N cm, 
and |M|=91.0Ncm. 



Problem 6.10: Forces in a three-leg stand 

Find the rod forces in a three-leg stand that is movably linked at the points A, B, C to a vertical wall 
and loaded at the point D by the force F. 

Solution Only longitudinal forces may act in the rods, because of the movable links of the suspension rods 
(neglect of bending forces). The forces at the cut-out branching point D are to be considered as 
external forces and are obtained from the equilibrium condition 

Fi + F 2 + F 3 + F = 0. (6.16) 

Using the unit vectors e, ( i = 1, 2, 3) along the rod axes and the magnitudes F, (i — 1 , 2, 3) of 
the rod forces, 6.16 may be written as 

Fiei T F 2 e 2 -1- F 3 e 3 = — F. (6.17) 

To get the rod forces, 6.17 is scalar-multiplied successively by the vectors e, x e, (i ^ j), where 
(e, x e.j) by definition points perpendicular to e,-, hence the scalar products e, ■ (e, x e ; ) vanish. 




The pair of forces responsible for the 
torque. 
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Using the definition of the triple scalar product A • (B x C), one then obtains from 6.17 for F t 
(i = 1,2,3) 

Fi = _ F.(e 2 xe 3 ) ; fi = _ F.(e i xe,) | ^ = _ F ■ (e, x e 2 ) _ (6 . 18) 

ei ■ (e 2 x e 3 ) e 2 • (e 3 x eO e 3 ■ (e 3 x e 2 ) 

Putting a coordinate frame into the branching point D according to the above figure, one gets for 
the unit vectors 



ei = (— cos a, sin a , 0) , 

e 2 = (cosa, sina, 0), ( 6 . 19 ) 

e 3 = (0, sin p, — cos P). 

Insertion of equation 6.19 into equation 6.18 yields 



F • (e 2 x e 3 ) = 



F x F y F z 
cosa sin a 0 

0 sin p — cos p 

— F x sin a cos p + F y cos a cos P + F z cosa sin p 



and 



ei • (e 2 x e 3 ) = 



— cos a sin a 
cosa sin a 

0 sin p 



0 

0 

— COS P 



= 2 sin a cosa cos p. 



From there one obtains for the component F\ 

F 1 Fy_ _ F z Unp \ 

2 \cosa sin a sin a ) 



( 6 . 20 ) 



( 6 . 21 ) 



( 6 . 22 ) 




z 




The forces acting at the three-leg 
stand. 
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The calculation of the scalar triple products for F 2 and F 3 from equation 6.18 runs in an analogous 
manner; one gets 



Fn= 1 (__F* Fy fitan/A 

2 \ cos a sin a sin a ) 



cos fi 



(6.23) 



Problem 6.1 1 : Total force and torque 

(a) Determine the components F y , F z of the force F = (2N, F y , F z ) acting at point Pi (1, 2, 3) m such 
that it is perpendicular to the plane defined by the three points Pi, P 2 {2, 3, 4) m, and P 3 ( 2, 2, 1) m. 

(b) What is the magnitude of the force F and which torque M does it apply with respect to the point 
P 4 (0, 1, 2) m? 

(cj What is the component of the torque vector M that points perpendicular to the plane? 




Solution (a) Because the vectors (rj — r 2 ) and (rj — r 3 ) are within the represented plane, the vector product 
R = (it — r 2 ) x (ri — r 3 ) yields a normal vector R perpendicular to the plane. If the force F shall be 
perpendicular to the plane, that is, parallel to the vector R, the following must hold: 

R = (R x , R y , R z ) = (n - r 2 ) x (r, - r 3 ) = XF = X (2 N, F y , F z ). (6.24) 

From there it follows that 

A=^tf r , Fy = { R y’ F * = \ R f (6-25) 

For the points i 3 ! (1,2,3) m, P 2 ( 2, 3, 4) m, P 3 ( 2, 2, 1) m, and / 4(C), 1, 2) m specified in the problem, 
one easily gets the position vectors r, (i = 1, 2, 3, 4) as well as their differences (r,- — r ; ) (i 7^ j): 



(ri - r 2 ) = (-1, -1, — 1) m, 
(n - r 3 ) = (-1, 0, 2) m, 

(n - r 4 ) = (1, 1 , 1) m. 



(6.26) 
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For the vector product in equation 6.24, it then results that 

R = (n - r 2 ) x (it - r 3 ) = (-2, 3, -l)nr = ( R x . R y , R z ). 

Inserting these values in 6.25, one obtains 

1 , , 

X = R x = -lm 2 r', 

2N 

F y = — = — 3N, (6.27) 

A 

R , 

F- = — = IN. 

Thus, the components of the force F are 

F = (2, —3. 1) N. (6.28) 

(b) The magnitude of the force F is obtained as 

|F| = (F 2 + F 2 + F 2 ) 1 ' 2 = (2 2 + 3 2 + 1) 1/2 N w 3.74 N. (6.29) 

The torque M with respect to the point P 4 results from the vector product 
M = (it — r 4 ) x F = (1, 1, 1) x (2, -3, l)Nm = (4, 1, -5) Nrn. (6.30) 

(c ) The component of the torque vector M perpendicular to the plane, that is, along the orientation of 
the force F, results from the definition of the triple scalar product as 

|M F | = [( ri -r 4 )xF]~=0. 



(6.31) 




7 Differentiation 
and Integration 
of Vectors 



Formation of the differential quotient: The vector A may occur as a function of a 

parameter. Let’s consider, for example, the position vector r(f) that — as a function of the 
time t — describes the path of a mass point. If one decomposes A into its components with 
respect to fixed unit vectors, then these components are functions of the parameter. We 
write 



A (u) — A x (u)e i + A v (n)e 2 + A z (u)e 3. 



(7.1) 



The differential quotient of a vector is formed by differentiating its components sepa- 
rately, as corresponds to the differentiation rule for sums. Because the unit vectors are not 
variables, they are conserved under differentiation, 

dA(u) A (m + Am) — A(u) 

= lim 

du Au—>o Au 

( A x (u + Au) — A x (u) A y (n + Au) — A v (m) 

= Inn — : ei + — - ; e 2 

Ah— >0 V All All 



+ 



A z (u + Au) — A z (u) 
Am 



-e 3 



The limit of the sum is equal to the sum of the limits, that is, when passing to the limit, 
one obtains 



dA(u) dA x {u) dAy(u) dA,(u ) 

— — = — e, + — ; e 2 + — e 3 . (7.2) 

au du du du 
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By comparing (7.1) with (7.2), one notices that the differentiation of a vector in an 
arbitrary coordinate frame with fixed unit vectors amounts to the differentiation of the 
components of the vector. Generally, the rule for the n-fold differentiation of a vector reads 



d n A(u) d n A x (u ) d n AJu) d n AJu ) 

= ei H : e2 H 1 — e3 

du n du n du n du" 



(7.3) 



Example 7.1 : Differentiation of a vector 

A (u) = (2 u 2 — 3zf)e! + (5 • cosz<)e 2 — (3 • sinz<)e 3 , 

A x (u ) Ay(u) A z (u ) 

= (2 u 2 — 3 u, 5 • cos u, —3 ■ sinzz), 

dA(u) 

= (4z< — 3)ei — (5 • sin ;z)e-> — (3 • cosw)e 3 

du 

= (4a — 3, —5 • sin u, —3 • cos u ), 

d 2 A(u) 

— — — = 4ei — 5 • cos u ■ e? + 3 • sin u ■ e ? 
du - 

= (4, —5 • cos u, 3 • sin u). 



For composite functions the usual differentiation rules apply. For example, for the product 
of a scalar function and a vector function, or for the scalar or vector product of two vector 
functions (parameter u), the product rule applies. 



Differentiation of the product of a scalar and a vector: 

d{(b(u)A(u)) d 

= —(<j>(u)A x (u)e i + (p(u)A v (u)e 2 + <j>(u)A z (u)e 3 ) • 

du du 

Now 

d(4>A x ) dtp . , dA x 

= — A x + <p 

du du du 

and analogously for the other components: 

-^(<M,) = (<P)Ai + -^-(Aj)(p 0 = 1, 2, 3). 

du du du 

This yields 



d((p(u) ■ A(w)) dtp dtp d<p dA x dA v dA , 

j = -t— A x e\ + — — A v e 2 + — A;.e 3 + (p — — ei + <p— —— e 3 + <p—^e 3 

du du du du du du du 



or simply 

d((p(u)A(u )) dp dA 

= — A + <p — ■ 

du du du 



(7.4) 
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Differentiation of the scalar product: One has 



d{A(u) ■ B(w)) d 



d 



du 



and therefore, 



= — ( y ^ Aj(u)Bj(u) = ^ —(Aj(u)Bj(u)) 



o'=i 



= E 



i=l 



dAj(u) 

du 



Bj(u) + A,(m) 



. , du 
dBi(u) 



du 



d(A(u) ■ B(m)) d A c/B 

= — B + A — . 

du du du 



(7.5) 



Differentiation of the vector product: It is performed analogously to the differentiation 

of the scalar product. Because the vector product is not commutative, one has to take care 
of the ordering of the factors. 



d 

— (A(n) x B(h)) 

tin 



dA(u ) 
du 



x B(n) + A(u) x 



dB(u) 

du 



(7.6) 



This is easily proved by checking the individual components (e.g., the x-component) on 
both sides of the equation. 



Example 7.2: Differentiation of the product of a scalar and a vector 

For the scalar function <p(x) = x + 5 and the vector function A(x) = (x 2 + 2x + 1, 2x, x + 2) the 
second derivative of the products ip • A is to be calculated. 

The differentiation of the product yields 



d 2 ((pA) d / dip dA\ 

dx 2 dx \dx ^ dx ) 



drip dip dA d 2 A 

dx 2 dx dx ^ dx 2 



The derivatives of the individual functions read 



dip d 2 ip 

dx ' dx 2 



dA 

dx 



= (2x + 2, 2. 1), 



d 2 A 
dx 2 



( 2 , 0 , 0 ). 



From the above, it results that 
d 2 ( ipA ) 

— I—L = (4x + 4, 4. 2) + (2x + 10, 0. 0) = (6x + 14, 4. 2). 
dx 2 



Application: Position, velocity, and acceleration of a mass point on a defined trajectory 

may be represented as vectors. The position vector for the motion of the mass point on 
an arbitrary trajectory B is the vector from the origin of the coordinate frame to the mass 
point; the variation of the position of the mass point with the time may be represented as 
time variation of the position vector (compare with the figure). 






42 



DIFFERENTIATION AND INTEGRATION OF VECTORS 7 




Definition of the orbital velocity: Ar = r(f + At) - r(t) is a secant vector to the orbit at the point r (f). The 
velocity is then given by v = lirriAf— ►o At /At. 



The velocity is defined as the first derivative of the position vector r(f) of the orbital 
curve with respect to the time: 



Ar r(f + At) — r(t) dr 

v = lim — = lim = — . 

Af^o A t Ar— >o At dt 



(7.7) 



From equation (7.7) one notices that the vector of the velocity represents the limit position 
of the secant through the position vectors r(t + At) and r (t) divided by the time interval 
At in the limit At 0, that is, the velocity points along the tangent to the trajectory at the 
point r (t). 

The acceleration is obtained as the first derivative of the velocity with respect to the time, 
or as the second derivative of the position vector with respect to the time: 



a (0 = 



d\(t ) 
dt 



= lim “ = 

Ar-s-0 At 



d(d r /dt) 
dt 



d 2 r{t) 

dt 2 



(7.8) 



Because the position vector is a vector, its derivatives with respect to the scalar time it) 
are again vectors. Thus, the velocity and acceleration are vectors, too. 



Problem 7.3: Velocity and acceleration on a space curve 

Let the position vector be given by r = ( t 3 + 2 1, —3e~', t ) m. Find the velocity and the acceleration 
as well as their magnitudes for the time points t = Os and t = 1 s. 

Solution For the velocity and acceleration, we get 

v(r) = f= (3f 2 + 2, 3e~‘ , 1) 

s 



a (f) = r = ( 6 1 , —3e ',0) — . 

s 2 

For the time t = 0, the results are 

v(0) = (2, 3, 1) — , a(0) = (0, -3,0)^, 

s s 2 
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/ 

v(0) = V14 

s 



m 

fl(0) = 3 — . 

s 2 



For t = 1 s, 



3 m 

v(l) = (5, 1) — , 

e s 



3 m 

a(l) = (6,— 

e s 2 



d(1) = 5.22 



m 



a( 1) = 6.1 



m 

T ' 



Example 7.4: Circular motion 

The Cartesian components of a circular motion are given by 

x(t) = R cos cut, 
y(t) = R sin cut, 
z(t) = 0. 



cu is the so-called angular velocity or also angular frequency. It is related to the revolution period T 
via u>T = 2n. The position vector is now 

r(0 = 0 ( 0 . y(t), z(t)) 

= Jf(t)ei + y(t)e 2 + z(t )e 3 , 

r(f) = (R ■ cos a>t, R ■ sin cut, 0) 

= R ■ cos a>t ei + R ■ sinnrf e 2 + 0e3. 

For the velocity one gets 
dr 

v = — = (— u>R ■ sin (ot, Ra> cos cut, 0) . 
dt 

There holds 



dr 

r ■ v = r • — =0 

dt 

=> v J_ r, 



for any time point, 




which is immediately clear for a circular orbit. 
For the magnitude of the velocity, one obtains 



u = |v| = ,/(— — ) + (-p") + ( _ r) — \/ cu 2 R 2 sin 2 cut + a> 2 R 2 cos 2 cut + 0 



dt 



dy 



dt 



dz 



dt 



/— ^ — 2n R circumference 

= v oi‘f? 2 (sin 2 cut + cos- cut) = cuR = = . 

T revolution period 

The acceleration is obtained as 
dy d 2 r 
dt dt 2 

= (—cu 2 R cos cut, —cu 2 R sin cut, 0) = —cu 2 (R cos cut, R sinaV, 0), 



= —cu r. 
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It turns out that the acceleration points opposite the orientation of the position vector (centripetal 
acceleration). The magnitude of the acceleration is given by 




Example 7.5: The motion on a helix 

The Cartesian coordinates of the helix read 

x(t) = Rcosoot, y(t) = R sincot, z(t) = bcot. 

The position vector is obtained by insert- 
ing in the relation 

r(f) = (x{t), y(t), z(t ) ) , 

that is, it holds that 

r(f) = (R cos cot, R sin cot, boot). 

Remark: b > 0 means right-handed he- 

lix, b < 0 means left-handed helix. 

The velocity results analogously to that 
of circular motion 

v = (— Rco sin cot, Rco cos cot, bco). 



1 1 Z 




The helix and its pitch. 



■> 

x 



The third component 1)3 = boo implies a uniform (constant) upward velocity (z-direction) if the 
parameter t represents the time. 

One has 

|v| = 1 / R 2 co 2 + b 2 oo 2 = io\J R 2 + b 2 , 

that is, the magnitude of the velocity is constant. 

The acceleration is the derivative of the velocity 

b = — or • (R cos cot, R sinoV, 0) = — co 2 r±, 
where 

r x = (R cos cot, R sin cot, 0) = (r • e r )e r 



and e, = (cos cot, sin cot, 0) is the polar unit vector in the x, y-plane. We thus obtain the same 
acceleration as for the circular motion. For the magnitude, it holds that |b| = oo 2 R. 



Integration of vectors: The integration rales may be applied also to vectors in the 

customary way. For a vector A depending on a parameter (e.g., u), it follows that 



J A (n) du 



I (A. v (n)ei + A y (u)e 2 + A z (u)e 3) du 
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= / A x (u)e\du+ / A v (m)C 2 <j?m+ / A z (u)e^du. 



If the unit vectors are constant, they may be pulled out before the integral symbol: 

J A(u)du = ei J A x (u)du + e 2 J A y (u)du+e 3 J A z (u)du. 

Thus, we may formulate the following rule: A vector is integrated by integrating its 

components. This vector integration graphically means a summation of a large number of 
vectors according to the integral limits; for example, the sum of all forces acting on a body. 
More strictly speaking: A (u) is a vector density, and d A — A(u) du is the vector associated 
with the interval du. These d A are summed to yield the integral. An example is the impulse 
of force K, which is understood as the force K acting on a body over a time interval; thus 
K = f At F(t')dt'. The impulse of force is therefore the sum of the forces V(t') acting 
during the time interval. For more details, see Chapter 17, equations (17.14) and (17.15). 



Example 7.6: Integration of a vector 



A = (2 u 2 — 3 11 , 5 cos u, —3 sin u), 



J A du = + CiJ e[ + (5sinu + C 2 )e 2 + (3cos« + C 3 )e 3 



(2 , 3 A 

= I -u — -u J e! + (5 sin«)e 2 + (3cosw)e 3 + Ciei + C 2 e 2 + C 3 e 3 



(2 , 3 2 \ 

= -u u ~ , 5 sin u , 3 cos u | + C. 

V 3 2 J 



The integration constants arising in the components are composed to the vector C. 



Problem 7.7: Integration of a vector 

Calculate 



A (n) dn with A = (3 re — 1, 2n — 3, 6 rf — An). 



Solution 



I A (n) dn = / (3 n~ — 1, 2n — 3, 6 n~ — An) dn 
Jo Jo 

= [( n 3 - n, n 2 - 3 n, 2 n 3 - 2 n 2 )]] 

= ( 6 , - 2 , 8 ). 
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Problem 7.8: Motion on a special space curve 

(a) Which curve is passed by the vector 

r(f) = ( x(t ), y(f), z(t)) = (f cos f, t sinf, f ) 
when f is running from 0 to 27r? 

(b) Calculate the velocity and acceleration of the point at the time f . 

(c) What are the velocity and acceleration for f = 0 and t = 2? 

(d) How do the magnitudes of radius vector, velocity, and acceleration vary for large time f? 



Solution (a) We first consider the vector r (f ) with l,(t) = 0 
(projection onto the x, y-plane). 

r (t) — ( f cos f , t sin f, 0). 

Because 

|r (f)| = (f 2 cos 2 1 + t 2 sin 2 t) 1/2 = t, 

there results a spiral line with a radius from 0 to 2?r. 

If z(t) — t additionally runs from 0 to 2n, we obtain 
a spiral line on the surface of a cone of height 2tc with 
the vortex at (0, 0, 0). 

The figure at the top of the facing page illustrates this 
result. 

(b) For the velocity v(t) and acceleration b (t), it re- 
sults that 

dr 

v(f) = — = (cos t — t sint, sin t + t cos t, 1), 
dt 

d\ d 2 r 

b W = FT = TT 
dt dt 2 

— (— sinf — sinf — f cosf, cosf + cosf — f sinf, 0), 
= (—2 sin f — f cos f, 2 cos f — t sinf, 0), 




The resulting spiral line with a radius 
varying from 0 to 2jt. 



(c) One has 




v(f = 0) = (1,0. 1); 


|v(f =0)|= a/5, 


v(f = 2) = (-2.23, 0.08, 1); 


|v(f=2)| = V6, 


b(f = 0) = (0, 2, 0); 


(N 

II 

O 

II 

»Q 


b(f = 2) = (-0.99, -2.65, 0); 


|oo 

II 

<N 

II 



(d) 



|r(f)| = (f 2 cos 2 f + f 2 sin 2 f + f 2 ) = \/2 |f | 



1/2 



|v(f)| = ((cos f — f sinf) 2 + (sinf + f cos f) 2 + l) 



I * 
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= (2 + t 2 ) 1/2 = |f| I 1 + — 



1/2 



= |fNl + i-0 



i»i 



|b(f)| = ((2 sin t + / cos t) 2 + (2 cos t 



= (4 + t 2 ) 1/2 = |f | ( l + -j 



1/2 



= |fNl + ^-0(l 



i»i 



\t\ (by a series expansion of the square root), 
t sin/) 2 ) 1 ^ 



|/| (by a series expansion of the square root). 



Problem 7.9: Airplane landing along a special space curve 

An airplane is landing. Thereby it is moving on the space curve 
r(f) = (x{t), y(t ), z(t)) = ( R cos cut, R sinmf, ( H — bcot)), 
with 



R = 1000 m, 




H = 400 m, 
b = H/6n, 
t e [0, 42 jt ] s. 
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Solution 



What is the velocity of the plane on landing (at t = 427r s)? Would you try landing this way? 



The velocity is calculated to be 

dr 

v = — = (— ( oR sintnf, coR cos cot, —bed), 
dt 

and its magnitude is 

|v| = ( co 2 R 2 sin 2 cot + co 2 R 2 cos 2 cot + b 2 co 2 ) 1/2 
= co(R 2 + b 2 ) 1 ' 2 . 



Obviously, it is independent of t ! Insertion of the values yields 



M = \ I looch + 



400 2 

(6tt) 2 



1/2 



sa 142.9 ms -1 = 514.4 km h _1 . 

This kind of landing is certainly unsuitable; the approach velocity should better be reduced. 




8 The Moving Trihedral 
(Accompanying 
Dreibein) — the 
Frenet Formulas 



In some cases it may be simpler to express velocity and acceleration in natural coordinates. 
This means that the velocity and acceleration are not derived from the variation of the 
position vector with the time, but from its variation with the passed way s, the arc length, 
the starting point being arbitrary. Let the curve itself be given by the position vector 
r(f) = (xi(f), X2(t), Xs(t)). For infinitesimally small segments, the increase of the arc 
length is |r/r| = ds. 




The magnitude of the velocity is 



M 



dr |r/r| ds 

dt dt dt 



In the limit of small At, the 
absolute value of the secant 
vector Ar becomes the line 
element ds, that is, |Ar| 
ds. 



In order to become independent of the coordinate frame, a set of orthogonal unit vectors 
is put at the point of the trajectory of the mass point given by s. The set of unit vectors 
moves along with the mass point: it is therefore also called the “moving trihedral” or 
“accompanying dreibein.” As unit vectors one uses 
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T tangent vector, 

N principal normal vector, 

B binormal vector. 

Because the vectors form an orthonormalized set, it holds that (N x B) = T, cyclically 
permutable. In the following we give the precise definition of these three base vectors of 
the moving trihedral and show how they are calculated for a given space curve r (f). 

The function r (?) describes a space curve depending on the time f as a parameter. To 
determine the moving trihedral, one has to convert the function r(f) into r(s); this is done 
by substituting the time t — t(s) from s = s(t) (compare with equation (8.1)). 

The moving trihedral is determined from the local properties of the trajectory, dr /ds is 
a vector along the limit position of the secant, i.e., the tangent. 




A curve in space and the moving trihedral (shown at an arbitrary point P of the curve). 



The magnitude of this vector is Ar|/ A.v. For infinitesimally small segments one has 
|z/r| = ds ; thus \dr\/ds = 1. Hence one has determined the tangent unit vector: 
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Because ds — \dr\, it also holds that 

dr dr / dt v 

|t/r| \dr / dt\ |v| 

In order to determine the principal normal 
vector, one first forms 

T T= 1. 



By differentiating the scalar product of the tan- 
gent vector, one obtains 



d 

ds 



d T 

(T • T) = T 

ds 



dT 

T =0. 

ds 




The difference vector AT is a measure of the 
curvature of the curve. AT points toward the 
“inner side” of the curve. A straight line has 
no curvature, and the normal vector N is not 
uniquely defined. 



Because the commutative law holds for the scalar product, T ■ d T /ds is zero. This implies 
that d T /ds is perpendicular to T. 

The vector c/T /ds gives the orientation of the principal normal vector. We characterize 
its position by constructing, besides the tangent defined by TO), a second tangent T(.v + A.y) 
(neighboring tangent) that differs from the first one only by an infinitesimal vector AT (see 
figure). The principal normal vector lies in the plane spanned by the two tangents TO) and 
TO + A,v). Because the magnitude of dT /ds in general differs from unity, one still has to 
introduce a factor k for normalization: 

m dT 

k • N = — , 
ds 

where k = \dT/ds\. This is the first Frenet formula. The value k is always defined as 
positive. This is possible since the orientation of N may be chosen in an appropriate 
manner. The factor k is called the curvature of the space curve. 

The third unit vector, the binormal vector, is formed out of T and N: 

B = T x N. 



The orientations of all three unit vectors are functions of the arc length. 

The vectors of the moving trihedral (accompanying dreibein) may be differentiated with 
respect to the arc length. The three differential quotients are called Frenet' s l formulas and 
read 



d T 

ds 



K N, 



d N 

ds 



rB — /cT, 



(8.2) 

(8.3) 



l Jean Frederic Frenet, b. Feb. 7, 1816 — d. June 12, 1900, Perigueux (Dordogne). In 1840 Frenet entered 
the Ecole Normale in Paris as a scholar, was appointed as professor in 1848, and taught until 1868 at Lyon 
University. His research mainly concerned problems of differential geometry, and in 1847 he found the Frenet 
forms independently of Serret. 
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= — rN or 

ds 



d B 

ds 



(8.4) 



r is a conversion factor and is called the torsion. The torsion describes the winding of the 
curve out of the T, N-plane. The quantity dB/ds is exactly a measure for this winding. 
From the curvature and torsion one gets 

1 1 

q — — curvature radius, a — — torsion radius. 
k x 



The curvature radius of a curve at a definite point equals the radius of the osculating 
circle having the same curvature as the curve at that point. 

Formula (8.2) has already been introduced as a definition. For a transparent derivation 
of the remaining formulas, one utilizes the statement that any vector may be represented as 
a linear combination of the three unit vectors. 



Derivation of the second Frenet formula: Because the moving trihedral spans the entire 

three-dimensional space, it holds that 

dN 

— = aT + /8N + y B , 
ds 

where a, /l, y are to be determined. Because N is a unit vector, N • N = 1. 

By differentiating the scalar product N • N, one obtains 

d d 

— (N • N) = — — (1) = 0 

ds ds 

or in other notation (product rule) d N / ds • N + N ■ d N /ds — 0, or, because the commutative 
law holds for the scalar product: 

dN 

2N = 0. 

ds 

Because in the nontrivial case neither N nor dN /ds is equal to zero, this means that 
d N /ds is perpendicular to N, i.e., there is no component of dN /ds along N. Therefore, 



dN 

= 0; i.e., — = aT + y B. 

ds 



(8.5) 



Moreover, according to the definition of the unit vectors, T • N = 0. By forming the first 
derivative of this scalar product one has 



d T dN 

N+T = 0. 

ds ds 



( 8 . 6 ) 



Using the first Frenet formula, we find 

dT 

N = kN ■ N = k. 

ds 



(8.7) 
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By inserting (8.7) in (8.6), we obtain 

dN dN 

k + T = 0, or T = -k. 

ds ds 

Multiplication of equation (8.5) by T yields 
dN 

T = aT • T + yB • T = a. 

ds r 

Because T • dN/ds = — k , it follows that a = —k. Hence 
dN 

— = — /cT + rB, 
ds 

where y = r is defined and inserted as conversion factor. 



Derivation of the third Frenet formula: We first try with the preceding trick and start 

from B • N = 0. If we differentiate the scalar product B • N, the product rule yields 

dB dN 

— -N + B- — =0. 
ds ds 

But this does not help immediately. We therefore simply start from the definition of B. 
Because B = T x N, it follows that 



d B d d T d N 

= (T x N) = — — x N + T x — . 
ds ds ds ds 

The first term of the equation may be transformed as follows: 

— xN = /rNxN = 0. 
ds 

By inserting equation (8.9) in (8.8), it follows that 

dB dN 

— = T x — . 

ds ds 

With 

dN 

— = rB — /rT, 
ds 

it follows that 

dB dB 

— = T x (rB — a:T), — = r(T x B) — k(T x T). 

ds ds 

Because 

TxB = — N and T x T = 0 , 



(8.8) 



( 8 . 9 ) 



it follows that 
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Darboux rotation vector: The Frenet formulas can be formulated in a very elegant way. 

To this end, we define a vector D as follows: 

D = rT + /cB. 

This vector D is called the Darboux 2 rotation vector. We now consider 

DxT = (rT + atB) x T 

= r(T x T) + r(B x T). 

Because T x T = 0 and B x T = N, it follows that 

D x T = kN. (8.10) 

Correspondingly, one has 

D x N = (rT + kB) x N 

= r(T x N) + k( B x N). 

Because B x N = — T and T x N = B, it follows that 

D x N = tB - n T (8.11) 

and there holds 

DxB = (rT + atB) x B 

= r(T x B) + at(B x B). 

Because BxB = 0 and T x B = — N, it follows that 



D x B = -rN. 



( 8 . 12 ) 



Using (8.10), (8.11), and (8.12), one may rewrite Frenet’s formulas in the following, 
highly symmetric form: 



d T 

ds 



= D x T, 



d N 

— = D x N, 
ds 



d B 

— = D x B. 
ds 



2 Jean Gaston Darboux , b. Aug. 14, 1842, Nimes — d. Feb. 23, 1917, Paris. Darboux came from modest 
relations. After graduating from Ecole Polytechnique and Ecole Normale in 1861, he decided for a teacher’s 
profession at the Ecole Normale. Supported by influential Parisian scientists, he got two teaching assignments 
after his doctorate in 1866. In 1881 he was appointed as professor. From 1880 on, he rendered merits as dean of 
the faculty of natural sciences on reorganizing the Sorbonne. From 1900 he served as permanent secretary of the 
Academie des Sciences. His main results concern the theory of areas. But he always aimed at joining to possibly 
all branches of mathematics, to penetrate them geometrically, and to work out the organic connection between 
mechanics, variational calculus, theory of partial differential equations, and theory of invariants. 
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Problem 8.1 : Curvature and torsion 

Prove the relation 

dr / d 2 r d 3 r\ r 

ds yrfi 2 ds 3 ) q 2 

Solution By inserting Frenet’s formulas and T = dr /ds, it follows that 



dr 


/ d 2 r d 3 r\ 




f dT 


/ rfN 


die y 


ds 


yrfi 2 ds 3 j 


= T 


X 
1 ^ 

1 


1 ^ 
^3 


+ dsV_ 



(dT dN\ 

= T • — x k — = T ■ (jcN x ic(rB - kT)) 

\ ds ds J 

= Tic 2 ■ (N x (rB - /cT)) = k 2 T • ((N x rB) - (N x acT)) 

= k 2 T • (N x rB) = k 2 rT ■ T = k 2 z = • 

Q 



Examples on Frenet’s formulas: 



Example 8.2: Frenet’s formulas for the circle 

Given the position vector 

r (f) = ( R cos cot, R sin cut, 0), 
calculate the vectors of the moving trihedral. 

Tangent vector: One has 

dr 

1 = ds' 

and with ds = |rfr|, it follows that 
dr dr/dt v 

\dr\ \dr/dt\ |v| 

The velocity is 

dr 




Tangent and normal vector of the circle. 



dt 

dr 

dt 



= v = (—Roosincot, Roocoscot, 0) = Rco{— sincuf, cos cot, 0), 



= |v| = cj R 2 a> 2 sin 2 cut + R 2 a> 2 cos 2 cot = Rco. 



Hence, for the tangent vector one obtains 
Ra>{— sin cur, cos ait, 0) 



T = 



= (— sin cot, cos cut, 0) . 



Rco 
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Normal vector: According to the first Frenet formula 

dT clT Idt 

k N = — = 

ds ds/dt 

We start from the time derivative: 

dT 

— = —co(coscot, sinaV, 0), 
dt 

ds dr 

— = — = Rco, 

dt dt 

dT — co(coscot, sinaV, 0) 

ds Rco 

This means 

1 r 

/cN = (coscot, sin cot, 0) = 

R R 2 

One has 

| (cos cut, sin cot, 0)| = 1. 

Because the curvature k is defined as a positive quantity, thus always k > 0, the following holds 

1 1 

*rN = — (— coscot, — smcot, 0), thus /c = |a:| |N| = — . 

R R 

Thus k = 1 / R, and consequently one has 
N = (— cos cot, — sin cut, 0). 

As was expected, the curvature radius q = \/k = R because R is the radius of the circle. For an 
arbitrary space curve the curvature radius in general varies continuously; it equals the radius of the 
osculating circle at a point of the curve. The geometric position of the centers of curvature of a curve 
(the centers of the osculating circles) is called the evolute. For the example of a circle the orientation 
of the normal vector is opposite to that of the position vector. The normal unit vector always points 
toward the center of the curvature circle. In this case the evolute is the center of the circle. 

Binormal vector: The vector B is calculated from B = T x N. 



B = — sin cot 



cos cot 0 = e 3 (sin 2 cot + cos 2 cot ) = e 3 . 



— cos cot — sin cot 0 



B = (0,0, 1). 



— = (0, 0, 0) = -rN, => t = 0. 
ds 

The torsion (winding) equals zero because the curve lies within a plane. One easily realizes that the 
torsion vanishes for all plane curves because T and N are within the plane, therefore B = T x N 1 to 
the plane and is therefore constant. Flence, from the third Frenet formula and from dB/ds — 0, it 
follows that r = 0. The torsion specifies how fast the curve is running out (winding out) of the plane. 
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Example 8.3: Moving trihedral and helix 

The moving trihedral of a helix is calculated analogously to the case of a circle. The position vector 
describing the helix in space reads 

r(f) = (R cos cot, R sin cut, beat). 



Tangent vector: 

^ v clr/dt (—Rco sin cot , Rco coscot, bco) 

|v| ds/dt J R 2 co 2 ( sin 2 cot + cos 2 cot) + b 2 co 2 

(-R sin cot, R cos cot, b) 

~ -JR 2 + b 2 ’ 

Normal vector: 

dT dT /dt — Rco (cos cot, sin cot, 0) 

— — = , ; = k|N. 

ds ds/dt *Jr 2 + b 2 ■ coj R 2 + b 2 

The curvature k is always defined as positive; correspondingly, the orientation of N is fixed 
(compare p. 51). One thus obtains 

R 

N = (— cos cot, — sin cot, 0), kl = — r -. 

R 2 + b 2 

The curvature of the helix is somewhat smaller than that of the circle, which is geometrically 
plausible. 



Binormal vector: One forms the cross product 

B = T x N. 

In determinant notation: 



B = 



1 



-Jr 2 + b 2 

b sin cot 



= ei 



-Jr 2 + b 2 

i 



Cl ^2 6.3 

— R sin cot R cos cot b 

— cos cot — sin cot 0 



—b cos cot (R sin 2 cot + R cos 2 cot) 

+ e 2 / „ + e 3 - 



-Jr 2 + b 2 



-Jr 2 + b 2 

(b sin cot, —b cos cut, R). 



-JR 2 + b 2 



For b — >• 0, B = constant = (0, 0, 1). To calculate the torsion, one forms 

c/B c/B /dt (1/V R 2 + b 2 )(bco cos cot, bco sin cot, 0) 
ds ds/dt cos/ R 2 + b 2 

b 

~ R 2 + b 2 



(cos cot, sin cot, 0), 
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d B 

ds 



-rN. 



The vector N has already been calculated above: N = (— cos cot , — sin cot, 0) . From there follows 
the torsion of the helix. It holds that 

—b b 

~ T = R 2 + b 2 ' T = R 2 + b 2 



The torsion radius: a = 1/r = (R 2 + b 2 )/b. 

For b = 0, it follows that r = 0. r is a measure for the variation of B. i.e., for d\\ /ds. In other 
words: r is a measure of how the curve is winding out of the plane. 



The three unit vectors T, N, and B define three planes that have particular names: 

T and N span the osculating plane, 

N and B span the normal plane, 

B and T span the rectifying plane. 



Remark: For a straight line r(t) = a+fe/c = 0 (q = oo) and r = 0 (a — oo). N and 

B may then be put arbitrarily _L to T = e. This is quite clear. 



Velocity and acceleration of a mass point on an arbitrary space curve: For arbitrary 

space curves it is sometimes convenient to express the velocity and the acceleration by 
means of the new unit vectors. After introducing the vector T one has 



T = — , 
M 



v = Ivl • T = vT. 



This relation may be used to derive the acceleration. 

d 2 r d\ d m dv m d T 

b = — — = — = — (vT) = — T + v—. 
dt 2 dt dt dt dt 

By transforming the second term, one obtains for the acceleration 



d T d T ds d T 

dt ds dt ds 



d r> 9 dT d i’ 
b=—T + v 2 — = zr T 
at ds dt 



9 dv v 2 
v 2 k N= — T + — N. 
dt o 



The acceleration is composed of two components: the tangential acceleration dv/dt T 
pointing in the tangential direction, and the centripetal acceleration v 2 /qN pointing toward 
the center of the circle of curvature. For a uniform motion of a mass point on a circle 
(Example 7.4) there exists only the centripetal acceleration, because dv/dt = 0 due to the 
uniformity of motion. 
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The E volute and the Evolvent: The e volute 

E it) of a curve r(r) is the geometric position of 
the centers of curvature (centers of the osculating 
circles) of the curve r(f): 

E(f) = r(r) + e(f)N(f) 

= r(t) + -!-n(o - 

K(t) 

where 1 /k = q is the curvature radius of the curve 
at point r. For plane cun’es, (r = 0) this holds: 
The tangents of the evolutes are simultaneously 
normals of the initial curve, because 




Parallel evolvents J(f) belonging to the 
curve r(f). 



d E 

dt 



dr 

dt 



57 N(,) + e 



d N 
dt 







ds 

dt 



(rB - kT) 



' ds 

= ti 57 



= ^N. 

dt 





(because r — 0) 



The evolvent (or involute or unwinding curve) }(t) is the geometric position of the arc 
length s plotted along the tangents: 

m = r (r) -s{t)- T(f). 



s is measured from an initial point Pq. Depending on the choice of the initial point Pq, 
one obtains a family of curves, whereby two evolvents in each point have a constant 
relative distance in normal direction. Such curves are called parallel curves. This is seen 
immediately by demonstrating that the tangent to the evolvent is perpendicular to the 
tangent of the initial curve, i.e., T • Ty = 0. But this is evident because 



dj ds ds m d T 

T j = — T T —s(t) — ; 

dt dt dt dt 



thus T 



j 



-N. 



=o 



If one is dealing with plane curves, then the construction of the evolutes and of the 
evolvents are in some kind of inverse relation with respect to each other. One finds: 



I. One of the evolvents of an evolute is the initial curve itself, symbolically written as 



Je f (s) = r CO- 
IL The evolute of each evolvent of a curve is the initial curve itself, that is. 



Ej r (0 = r CO- 
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Here we have written the corresponding initial curve as an index; thus J r (.y) is the 
evolvent of the curve r(.y), and Ej r (,y) is again the evolute of the evolvent J r (s). 

We prove the second assertion: It reads 

E Jr (s) = Jr (.S') + — Nj = (r(.s) - >sT r (s)) + — Nj . 

Kj Kj 



The normal of the curve J r is obtained by differentiation of the tangent vector Tj = — N r 
with respect to the arc length s j (i.e. not with respect to .v = s , !). Therefore, 



Nj 



l^Tj _ 

Kj dsj 
1 ds 
kj dsj 



1 ds /-</N r \ 

Kj dsj V ds ) 

(TB r -KT r )= — -^Tr 

Kj dSj 



if the torsion r vanishes (plane curve!). 

The derivative of the arc length of the curve r with respect to the arc length of the evolvent 
J r is obtained because 



d Jr 
dsj 



d J r ds 
ds dsj 1 



thus 



ds 

dsj 



|Tj| _ 1 _ 

\dj r /ds\ s\dT r /ds\ sk 



Because Nj and T r are unit vectors, it must hold that 



k ds ds 11 

— = 1 or Kj — K- — —K — = 

Kj dSj dSj SK S 



We see that the curvature radius of the evolvent just equals the corresponding arc length 
s of the “unwound” curve, as is expected clearly. 

For the evolute of the evolvent, we now obtain 

Ej r = r(.S’) - sTr(j) + — Nj 
kj 

= r(.s’) - jT r (s) + — T r (s) 
kj 

— r(.y) - sT r (s) + sT r (i) = r(i). 



Thus, assertion II is proved. Assertion I may be proved in a similar way. 



Note: By adding a term pointing along the binormal direction, the definition of the evolute 

may be generalized in such a way that the assertion holds also for general space curves 
with torsion r ^ 0 (compare to Example 8.6). 




EXAMPLES ON FRENET'S FORMULAS 



61 



Example 8.4: Evolvent of a circle 

The evolvent of a circle is a spiral. The centers of curvature of this spiral are located on the circle, 
which therefore is the evolute of the spiral (compare the figure). 




Problem 8.5: Arc length 

Calculate the arc length of the space curve given by 
r(f) = 3 cosh(20 e, + 3 sinh(20 e v + 6fe, 
for the interval 0 < t < tc . Outline the curve! 



Solution 



One has 

i = 

dr 

dt 

dr 

dt 



f 


f ds 


■ 1 dr 


r ,= j 


ht dt = ) 


\ dt 



because ds = |dr|, 

6 sinh(2r) e v + 6 cosh(2r) e-y + 6e;., 

6 yj sinh 2 (20 + cosh 2 (2f) + 1 = 6-J 2 cosh 2 (2f), 




because sinh 2 jf = coslr x — 1 — >• \dr /dt \ = 6V2cosh(2r), 

7T 2 JT 

s = J 6\/2cosh(2 1 ) dt = -6V2 J cosh.rdx = 3\/2 sinh(2^). 



The space curve comes from the first octant, intersects the x, y-plane at the point (3,0,0), and enters 
the eighth octant — twisted hyperbola: Consider the x , y-components: x = cosh 2 1, y = sinh 2 1, then 
form x 2 — y 2 = cosh 2 2 1 — sinh 2 2t = 1 — >■ x 2 — y 2 = k ( k = constant). This is the equation of a 
hyperbola, see Chapter 26. 
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Example 8.6: Generalization of the Evolute 

The definition of the evolute may be extended to the case of nonplanar curves, i.e., curves with torsion 
r(s) ^ 0 in such a way that 

J,, r Cv) = r(s) (8.13) 

further holds. For this purpose we start from a general ansatz allowing that the evolute runs out of 
the osculating plane of the curve r(s), namely 

E(s) = r(,s) + A(,s)N(s) + /x(*)B(s) (8.14) 



with two indeterminate functions A(,y) and /x(s). 

To calculate the evolvent of E, one needs the derivative 



d E 

ds 



dr dX d^i da dB 

= 1 N + A 1 B + /x — 

ds ds ds ds ds 



= T(1 — kX) + N(A — fir) + B(/x + tA), 



(8.15) 



where the Frenet formulas have been utilized. The dot denotes differentiation with respect to s. 
The evolvent of the evolute then has the form 



JetCs) = E r (s) - 5 £ (s)T e (s) = E r (s) - 

ds E 

ds dE 

= r + AN + /xB — s E — — 

CIS £ CIS 



= r — j £ T(1 - kX) + N 

ds E 



ds .. . 

A - s E ~ — (A - /xr) 
ds E 



+ B 




ds 

Se — (m + tA) 
ds E 



(8.16) 



In order to fulfill 8.13, all additional terms on the right side of 8.16 must vanish. Because the 
vectors of the moving trihedral are orthogonal, one is led to three independent equations: 



1 - kX = 0; 

ds ■ 

A — s E — — (A — fir) = 0; 
ds E 

ds 

/x — s E (/x + rA) = 0. 

ds E 

The first equation again yields the old result 



A(s) = 



1 

k(s)' 



We now resolve equation 8.18: 



ds X 

ds E X — /xt 



1/k _ -k 

( — 1 //C 2 )/F — /XT k + /XT/C 2 



(8.17) 

(8.18) 

(8.19) 



( 8 . 20 ) 



( 8 . 21 ) 
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and insert this in equation 8.19: 

M + v— * j (A + rk) = 0. 

K + flTK ~ 

This is a differential equation of first order for the function fi(s), 
, k t 

A + r K[i H — /t H — = 0. 

K K 

In order to solve 8.23, we multiply by k, 

( /cA + + r/c 2 /t 2 T t = 0. 

We substitute T(i) = k(s)/j.(s); hence 

— Y + t(Y 2 + 1) = 0. 
ds 

This may be integrated by separation of variables. 



7t^T- + / — 



+arccotT = 



S 

j ds'x{s') 



fj.(s ) = cot f f ds' r(s') + C I . 

" (5) 1 / ) 

The generalized definition of the evolute therefore reads 



1 1 

EO) = r(s) + — -N(s) + — — cot 
k(s) k(s) 



s 

J ds' t(s') 



) + C B(s) . 



Because C is an arbitrary constant, there exists an entire set of evolutes. 




Q Surfaces in 



It may happen that the position vector is not a func- 
tion of one parameter only but depends on two pa- 
rameters u and v: 

r ( m, v) = (x(u, v), y(u , v), z(u, v)). 

The position vector then describes a surface in 
space. This shall be visualized: Let r be a func- 
tion of two parameters u and v. We first choose a 
fixed value i>i for v and let u vary continuously, 
r (u, Ui) then describes a space curve (compare to 
the figure). 

Now we choose another fixed value of v that is 
not widely spaced from V\ and denote it by ih- u is 
again varied continuously. 

There results a space curve r(u, ih) that does not differ too much from r(u. ui). This 
procedure may be repeated many times, and one obtains many neighboring space curves 
(see figure overleaf). 

Then, the same procedure may be performed in the opposite manner. By choosing a fixed 
value for u and varying v continuously, one obtains distinct neighboring lines r (u„, v) for 
a fixed u n (see next figure). 

If the spacings between u and v become more and more dense, one obtains a surface 
in space. One may form the derivative along such a curve (e.g., fixed u — h 2 and varying 
v). The derivative in which one of the parameters is considered variable while the other 
parameters are considered fixed is called a partial derivative and is denoted by a round 3 
(spoken: "d partial” or “d partially derived with respect to”). 

dr(uj,v) dr (u,v) 

u — ui = constant : = r„ = . 

dv dv 



Space 




Illustration of the space curves r(u, v n ). 
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In the same way one forms the tangent vector r u : 

dr(u, Vi) 3r (u, v) 

v — Vi = constant : = r u = . 

du 3 u 

The plane fixed by r M and r„ is called the tangent plane of the surface. From r„ and r„, one 
easily forms the normal vector n, which is perpendicular to the tangent plane. 
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The normal vector is 



n(«, v) = 



r„ x r„ 
|r„ x r„ 




If r„ • r„ = 0 at any point of the surface, the mesh 
formed by the curves for u = constant and v = 
constant, respectively, is called an orthogonal mesh. 

For example, the meridians and parallels of constant 
latitude of a sphere form an orthogonal mesh. A sur- 
face for which a normal vector may be constructed 
at any point is called orientable. There are surfaces 
with only one side, as for example the Mobius strip 
(see Section 14). On such a surface, any point can be 
reached from any other point by a continuous dis- 
placement of the normal vector. Such surfaces are 

not orientable. Orientable surfaces have inner and outer sides. By a continuous displacement 
of the normal vector, one stays always on the same side of a orientable surface. The normal 
n of an orientable surface is defined as positive for external (convex) surfaces, and negative 
for concave ones. 



Illustration of a surface in space with 
tangent and normal vectors and a tan- 
gent plane in one point of the surface. 



Example 9.1 : Normal vector of a surface in space 

The position vector r(u, v) = a cos u sin v ei + a sin u sin v £2 + a cos v e 3 with variable parameters 
describes a surface in space. 

Find the normal vector as a function of u and v. 



Solution r„ = —a sin u sin v 6| + a cos u sin v e? + 0 e? . 

r„ = a cos u cos v ej + a sin u cos v ei — a sin v , 



ei 62 63 



r„ x r„ = 



|r„ x r„[ = 



n = 



—a sin u sin v a cos u sin v 0 
a cos u cos v a sin u cos v —a sin v 

—a 2 cos u sin 2 v — a 2 sin u sin 2 v — a 2 sin v cos v , 

a 2 \fcos 2 u sin 4 v + sin 2 u sin 4 v + sin 2 v cos 2 v 
a 2 ^/lcos 2 u + sin 2 u ) sin 4 v + sin 2 v cos 2 v 
a 2 J sin 2 v (sin 2 v + cos 2 v) 
a 2 \ sin u|, 

(— cos u sin v, — sin u sin v, — cos v) for sin v > 0. 
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The result means that the normal vector always points opposite to the position vector, which is the 

case for a sphere. One may easily prove that the function of the position vector represents a sphere, 

by calculating the magnitude of the position vector: 

x = a cos u sin v, 
y = a sin u sin v, 
z = a cos v . 

The absolute value (the normalization) |r[ = r of the position vector is calculated from 
|r | 2 = x 2 + y 2 + z 2 = a 2 ( cos 2 u sin 2 v + sin 2 u sin 2 v + cos 2 u) 

= a~(sin t> (cos u + sin‘ u) + cos v) 

= a 2 . 

From there it follows that r = a — constant, i.e., the given position vector determines the surface 
of a sphere. 

Because 

r„ • r„ = —a 2 sin u cos u sin v cos v + a 2 sin u cos u sin v cos v + 0 (—a sin v ) = 0, 

the mesh spanned by the u, u-lines represents orthogonal coordinates. One easily confirms that the 
n-u-lines are the meridians and parallels of equal latitude on a sphere. 




- 4 /\ C oor d' nate 

I w Frames 



In an n -dimensional space one may always define n linearly independent base vectors out 
of which any arbitrary vector may be composed by a linear combination. For the sake of 
simplicity, vectors of magnitude unity are usually adopted as base vectors. 

Corresponding to the number of base vectors, the position of an arbitrary point may 
be specified by n independent real numbers i = 1 Each coordinate frame 

is characterized by a mutually unique assignment between the space points and these n 
numbers, the coordinates. 

A vector in the n -dimensional space reads 

n 

r = 7>e ; , 

! = 1 

where the n base vectors e, again shall satisfy the orthonormality relation e, ■ e ; = <5, j. 
The scalar product of two n -dimensional vectors a = {a, } and b = { b/ } may be defined by 
a • b = 1 a ibj, in analogy to the three-dimensional space. 

The introduction of a coordinate frame implies that the coordinates of a space-fixed point 
change if the frame is displaced or rotated. From there it follows that for any special system 
a reference point and a definite orientation in space must be given. 

Physically seen, both quantities may be fixed by tying the coordinate frame, for example, 
in a rigid body as a reference body; in a completely empty space it would make no sense to 
speak of the position of a point. Of course, a coordinate frame must not be “at rest” (e.g., 
all frames tied to the earth are accelerated frames due to the earth’s rotation). 

Special examples 

1. The position of a point on an arbitrarily curved line (n = 1) 
is already specified by giving one number. In the simplest case 
one adopts the arc length s measured from a reference point in 
a defined direction of motion as a “natural parameter.” This is a 
one-dimensional space. 




A caterpillar on a blade 
of grass. 
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2. The surface of the earth, although being formed in 
a highly complicated manner (mountains, etc.), is an 
area with n = 2. Each point on it may thus be uniquely 
determined by two numbers. As is known this may be 
achieved by fixing two angular quantities: geogra- 

phic length and latitude. Arbitrarily chosen reference 
quantities are the zero meridian through Greenwich 
(geographic length = 0) and the equator (geographic 
latitude = 0). This is a two-dimensional space. 

In order to change from one coordinate frame (q \ , qi , q 3 ) to another one (here specifically 
the Cartesian frame: x, y, z), the following equations have to be set up: 

Transformation equations: 

qi=qi(x,y,z) x = x(qi, q 2 , q 3 ) 

qi = < 72 U'- y, Z ) and their inversion y — y(q \ , c/ 2 , q 3 ) ( 10 . 1 ) 

q 3 =q 3 (x,y,z ) Z = z(qi, qi, qi). 

Cartesian coordinates: Given are the three base 

vectors e1.e2.e3 along the directions of three mu- 
tually perpendicular axes. The coordinates x, y, z 
of a point P are the projections of the position z 

vector r — OP onto the axes, 

r = xei + ye 2 + ze 3 , |e,j = 1. 

By convention the three unit vectors form a 

right-handed frame. Because they are mutually 

perpendicular, they constitute an orthogonal frame. 

Moreover, the unit vectors are always parallel to 

the axes, that is, fully independent of the position 

,, . ,, . The definition of Cartesian coordinates, 

of the point P in space. 

This constancy of direction of the unit vectors combined with their orthogonality is 
the reason for preferred usage of Cartesian coordinates. For many special problems with 
particular symmetry, it turns out as convenient to use coordinate frames that are adapted to 
the geometric conditions and therefore simplify the calculations. For example, the motion 
of a plane pendulum may be described in terms of one angular coordinate, the motion of a 
spherical pendulum in terms of two angular quantities. 

Curvilinear coordinate frames: To explain the denotation, we suppose the coordinates 

x, y, z of r to be expressed by q\, q 3 , q 3 according to (10.1). There results 





r(<7t, <72, < 73 ) = [x(qi,q2,q3),y(qi,q2,q3),z(quq2,q3)} ■ 
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Two of these three coordinates q\, q 2 , q 3 are now kept constant; let only the third one be 
variable. All points satisfying this condition are located on a space curve. There arise the 
three coordinate lines: 

Li : r = r(qi , q 2 = c 2 , q 3 — c 3 ), 

L 2 \ r = T{q x = ci,q 2 ,q 3 =c 3 ), (10.2) 

L 3 : r = r(qx = ci, q 2 = c 2 , q 2 ). 

As is immediately seen from the scheme, the three coordinate lines have exactly one 
common intersection point P (c\, c 2 , c 2 ). 

In the Cartesian frame these lines are straight lines parallel to the three axes. If, however, 
at least one of the lines is not straight, one speaks of curvilinear coordinates. One may still 
proceed one step further and keep only one of the three coordinates constant, while the 
other two remain variable. There arise two-dimensional (in general curved) areas in space. 

Coordinate areas: 



Fi : 


r = r(qi 


— Ci, q 2 


,q 3 ), 




F 2 : 


r = r( qi 


,q 2 = 


c 2 ,q 3 ), 


(10.3) 


F 3 : 


r = r(qi 


,q 2 


, c/3 = c 3 ). 





One may imagine the coordinate lines 
as resulting from the intersection of two of 
these areas. In the Cartesian system the co- 
ordinate areas are planes with the common 
point P . 

Generally an arbitrary point may be rep- 
resented as the intersection point of its three 
coordinate areas (and, of course, also co- 
ordinate lines). One presupposes that each 
space point is traversed by exactly one area 
from each of the three sets of coordinate 
areas. The three fixed parameters of these 
areas are the coordinates of the point. 

The vector r(r/i , q 2 , q 3 ) as a function of the three parameters q\,q 2 , q 3 describes a space 
region. Actually, if one of the coordinates is kept fixed, e.g., q 3 = q 3 , according to Chapter 9 
we are dealing with an area in space. If q 3 changes to q 3 = q 3 + A q 3 , a neighboring area 
emerges. If q 3 is running continuously, there emerge more and more arbitrarily densely 
located areas in space that, in total, cover a space region. 




Illustration of coordinate surfaces. 



General specification of base vectors: As normalized base vector (unit vector) e 9l at 

point P , we choose a vector of magnitude 1 tangential to the coordinate line q 2 = c 2 , q 3 = c 3 
at P . Its orientation shall correspond to the direction of passage of the coordinate line with 
increasing value q\. 
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This introduction of the unit vector corresponds exactly to the geometric meaning of 
the partial derivative; hence, e 9l may be calculated by partial differentiation of the position 
vector with respect to q\ and subsequent normalization: 

dr /dq\ 3r 3r 

e gl = ;a , or —=h l e qi or — = A,-e 9i ; i = 1,2,3. (10.4) 

|3r/3<?i| dq\ dqi 

Here h, are scaling factors, namely hi = 3r/ik/, |. In curvilinear coordinate frames the 
direction of at least one of the coordinate lines changes by definition. Therefore, these 
frames are, contrary to the Cartesian frame, coordinate frames with variable unit vectors. 



Cylinder coordinates: The coor- 

dinates used are 

< p : angle between the projection of 
the position vector onto the x, y- 
plane and the x-axis, 

q : separation of the point from the 
z-axis, 

z: length of the projection of the po- 
sition vector onto the z-axis (as in 
the Cartesian frame). 

The coordinate areas extend to in- 
finity (see figure, showing limited 
sections) and are 

q — Q \ : circular cylinders about the z-axis, 

<p = cpx : half-planes containing the z-axis, (1 0.5) 

Z = Zi : planes parallel to the x, y-planc. 

Coordinate lines are two straight lines and a circle. 




The definition of cylindrical coordinates. 



Transformation equations: From the figure one may directly read off the relations: 

r = (xi, JC 2 , X 3 ) = ( q cos <p, q sin <p, z) 



or in detail: 

x — q cos (p , q = yj x~ T y 2 , 

v y 

y — q sin</), (p — arctan — = arcsin — , (10.6) 

X Q 



z = z, 



z = z. 
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To ensure that one point cannot be characterized by distinct combinations of coordinates, 
we agree on the following restrictions: 

q > 0; 0 < <p < 2 7t. 

The representation is not completely unique since the angle remains indefinite for points 
with q = 0. But inversely — and this is the more important requirement — to each triple 
p, </?, z only one space point is associated. 



Unit vectors: According to the geometric 

introduction as tangent vectors to the coordi- 
nate lines, the unit vectors e e , e r are given 
by 



3r/3 q (cosy?, sin q>, 0) 



— 



e 7 = 



|3r/3p| 
3r /3 (p 
|3r/3^| 
3r/3 z 
|3r /dz\ 



1 

(— siny), cos <p, 0) 



0 



(0,0.1) 

1 



( 10 . 7 ) 



e e is parallel to the x, y-plane and points in 
radial direction from the "-axis outward. 



is tangent to the circle z = Zi, Q = Q i, 
that is, also parallel to the x, y-plane. 



Z 




X 

Illustration of cylindrical coordinates. 



e, corresponds to the Cartesian e?. 



Thus, e e and e lf> may be projected onto the x, y-plane without any changes. One has 
e e = cosy?ei + simpe 2 . 



= cos y<p + — j ei + sin [cp + — j e 2 = - sin^ei + cosy>e 2 , 
e e = (cos (p, sin <p, 0), 

= (— simp, cosy?, 0), ( 10 . 8 ) 

e- = (0, 0, 1). 

The same result follows by partial differentiation of r with respect to Q,<p,z and subsequent 
normalization (see equation (10.7)). 

To check the unit vectors, we form the triple scalar product 



cos (p siny? 0 

— sin <p cos <p 0 
0 0 1 



C- ) — 



= 1. 
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This is the unit volume spanned by the vectors e g , e v , e ; . Thus, the cylindrical coordinates 
form an orthogonal frame with variable unit vectors. For solving kinematic problems it is 
important to know the derivative of the unit vectors with respect to time. Let the functions 
g(t), (pit), z(t ) be known. The generalization of the chain rule for a function of several 
variables then yields 



de e 3e e dg de e dtp de e dz 

dt dg dt dtp dt dz dt 

— 0 + (— sin^>, coscp, 0)tp + 0 =tpe ip , 



i 

— - = (— cos tp, — sin (p, 0 )tp 
dt 



~(pe„ 



da, 

= 0 . 
dt 



( 10 . 9 ) 



The derivative of a vector e of constant magnitude has no component along the direction 
of e and hence must be perpendicular to it: e • e = constant => e • da/ dt — 0! 

The equations given above fulfill this condition! We still note that from now on we shall 
frequently abbreviate the time derivative of a quantity by a dot above this quantity, as, for 
example, dtp/dt = tp or de B /dt = e e , etc. 



Velocity and acceleration in cylindrical coordinates: Let a point move along a path 

described by the position vector r (?). One has 

(a) the velocity v(r) = dr/dt, 

(b) the acceleration b(t) = d 2 r/dt 2 = d\/dt. 

In cylindrical coordinates let g(t), tp(t), z(t) be given. The position vector is 

r = ge e + ze z . (10.10) 

Note: These base vectors are now not fixed but are coordinate-dependent by themselves. 

One has to take care in component representation: For instance one cannot simply dif- 
ferentiate r = (g, 0, z) ! In order to avoid errors, one has to write out the vector, as, for 
example, 

(a) 

r = ge e + ge e + ze z + ze z . 

This yields the velocity. 
r = ge e + g<pa v + ze z . 

(b) 

r = (ge e + ge e ) + (gtpe^ + gipa^ + gtpa^) + (lie- + ze z ). 



( 10 . 11 ) 
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This yields the acceleration: 

r = (e — Q<p 2 )e s + {.Qtp + 2 Q(p)^ + ze ; . (10.12) 

Hence, in the cylindric frame both the velocity and acceleration are composed of three 
components: a radial component, an azimuthal component, and a component in the z- 

direction. 



Spherical coordinates According to the figure below, the coordinates are 



r : length of the position vector, 

f) : angle between the position vector and the z-axis (polar angle), 



cp: azimuth (as in the cylindric frame). 



The previous figure illustrates the various coor- 
dinate areas and lines. The point P is the inter- 
section point of a circular cone about the z-axis 
with the vortex at the origin 0, a half-plane includ- 
ing the z-axis, and a sphere with the center at 0 
that results by keeping the radius r constant and 
varying the two angles. The coordinate lines are 
two circles and a straight line: (1) r = constant, 
q> = constant, i) variable — meridian; (2) r = 
constant, ft = constant, <p variable — parallel of 
constant latitude; ( 3 ) <p — constant, ft = constant, 
r variable — radial ray. 

The coordinate areas are a conical area (ft = 
constant), a half-plane (<p — constant), and a 
spherical area (r = constant). 



Z 




The definition of spherical coordinates. 



Transformation equations 

r = jciei + X2ez + -*363 = r sin ■& cos cpe \ + r sin ft sirups + r cos ft e3. 

When the equations are written in detail, we get 

x — r sin ft cos cp, r = ^ x 1 + y 2 + z 2 , 

y 

y = r sin ft sm<p, ^ = arctan — , (10.13) 

x 

ft — arctan x 1 + y 2 /zj ■ 



Z = r cos ft, 
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Coordinate surfaces and coordinate lines for cylindrical coordinates (left) and spherical coordinates (right). 
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To reach uniqueness, the following restrictions are agreed upon: 
r > 0, 0 < y < 2i r, 0 < # < n. 



( 10 . 14 ) 



Unit vectors for spherical coordinates: The position vector is 



r = r(sin# cos y, sin$ sin<p, cos $). 



Partial differentiation yields 



— = (sin cos y, sin i? sin y, cos #), 
dr 



— = r(cos f) cos y , cos 1 } sin^), — sin ft), 
dtf 



— = r(— sin & sin y, sin ft cos y, 0), 
d y 

The unit vectors follow by normalization: 

e r = (sin & cos y, sin i) simp, cos ft); 
e§ — (cos ft cos y, cos & sin y, — sin i})\ (10.16) 

e v = (— simp, cosy, 0). 



h r — 



h» — 



h,„ — 



3r 


= 1, 


dr 


dr 




d& 


= r. 


dr 


— r sin #. 


dy 



( 10 . 15 ) 



Geometrical interpretation: One has re, = r; hence, e, points along the position vector, 

that is, it is the normal to the surface of the sphere. 

lies tangential to the circle r = ri, ft = &i, namely, parallel to the x, y-plane. Its 
component representation may accordingly be seen from the previous figure when setting 
the circle radius equal to r sin i) . 

e„ has a component sin i) along the negative "-direction. We know that e„ is the tangent 
vector of the ^-coordinate line, namely the tangent to the meridian. The question whether e,, 
points upward or downward is decided by the /-component (— sin ft): e„ points downward 
as in the figure. 

One may easily convince oneself that the spherical coordinates also constitute an 
orthogonal frame with variable unit vectors, by evaluating the triple scalar product 
e,. • (e# x e v ) = 1. 

We write the equations (10.16) explicitly: 

e ( . = sin ft cos y ei + sin ft sin ye 2 + cos & e 3 ; 

e# — cos & cos ye 1 + cos & sin ye 2 — sin ft e 3 ; 

e v = — simpei + cos<pe 2 + 0 e 3 . 



( 10 . 17 ) 
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and solve them for ei, e 2 , e? according to Cramer’s rule 1 . For example, for ei one finds 





C r 


sin d sin <p 


cos d 


/ 


sin d cos <p 


sin d sin <p 


cos d 


ei = 


e# 


cos d sin <p 


— sin d 


/ 


cos d cos (p 


cos d sin <p 


— sin d 




£(p 


cos <p 


0 


/ 


— sin^ 


cos q> 


0 



e,- sin d cos ip + e$ cos d cos cp + (— sin tp) 

sin 2 d cos 2 q> + cos 2 d cos 2 tp + sin 2 (p 
ei = sin d cos q> e,. + cos ft cos q> e# — sin <p , 

and similarly for e 2 and e 3 : 

e 2 = sin d sin <p e, + cos d sin tpe# + cos <p e p ; 
e 3 = cos d e,- — sin d e,-> . 



(10.18) 



(10.19) 



Velocity and acceleration in spherical coordinates: To calculate the velocity and ac- 

celeration in spherical coordinates, we still need the time derivatives e, , e,> . e (/ , . One finds 



3e,- • 3e,- . 

e,. = — H <p 

3d 3 y 

= (cos d cos (p , cos d sin <p, 
= de# + sind^e^, 



■ sin d)d + (— sin d sin<p, sin d cos cp, 0 )<p 



( 10 . 20 ) 



and similarly 

e f) — — d e r + cos d (pZy, 

= — sin ft <p e,- — cos d <p e#. (10.21) 

Now we may calculate the velocity in spherical coordinates. The following hold: 



r = re,, 
r = re,. + re r 

= re,- + rde# + r sin ft £>e p , ( 10 . 22 ) 

r = re,. + re r + rdep + rdep + rde# 

+ r sin ?? e y + r cos d d <p + r sin d ^3 + r sin d <p 

(after inserting (10.20) and (10.21) ) 



1 Gabriel Cramer , b. July 31, 1704, Geneva, as son of a physician — d. Jan. 4. 1752, Bagnols near Ntmes. After 
his studies at the university of Geneva, Cramer became appointed as professor for philosophy and mathematics. 
From 1727-1729 he made an informative trip through many European countries. After his return home Cramer 
held important municipal posts in Geneva. His rapidly decaying state of health led him to southern France where 
he soon died. His main work is the Introduction a 1’ Analyse des Lignes Courbes Algebriques (1750), where among 
other things the theory of solving systems of equations by means of determinants is outlined. 
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b(p 



= b r e r + + b^e^. ( 10 . 23 ) 

If # = 7r/2, that is, sin# = 1, # — 0, cos # = 0, (10.22) and (10.23) turn into 
r = re r + rip^ 
and 

r — (r - r<p 2 )e r + ( 2 rip + r<p)e v , 

respectively. These expressions for velocity and acceleration in plane polar coordinates are 
already known from the discussion on cylinder coordinates. 



Problem 10.1: Velocity and acceleration in cylindrical coordinates 

A particle moves with constant velocity v along the heart curve 
or cardioid r = k(l + cosp) (Greek kardia = heart). Find the 
acceleration a, its magnitude, and the angular velocity. (Note 
that r denotes here the coordinate g of the cylindrical coordinate 
frame.) 



cos ip = 



1 and sin <p = 1 — ( 



-•-(i-'V- 



2 r 
= 2 - 
k 



r~ 

¥ 



According to 10.25 we obtain 
r 2 = k 2 sin 2 (pip 2 = k 2 



r r 
2 k~ ¥ 



( p 2 = 2 krtjr — r 2 (p 2 



Solution The differentiation of the path equation with respect to time yields 

r = k{\ + cos <p), ( 10 . 24 ) 

r = —k sin <p <p, ( 10 . 25 ) 

f = — k((p 2 cos <p + ip simp). ( 10 . 26 ) 

For the discussion below it is useful to conclude from 10.24 that 

,2 




( 10 . 27 ) 



( 10 . 28 ) 



Because we are dealing with plane polar coordinates, we write for the radius vector 
r = re, , 



r = re r + r<pe v , 

r = (r — r<p 2 )t r + (rip + 2 r<p)e v . 



( 10 . 29 ) 

( 10 . 30 ) 

( 10 . 31 ) 
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Because the velocity is given as constant, from 10.30 it follows that 
v = y/ r 2 + r 2 <p 2 , 

and with 10.28 it follows for the angular velocity that 

7 Wr' 

because namely 

u= (/( 2 ^“S)* 2+r2 )^ = ^ : 

hence 

<f> = " • ( 10 . 32 ) 

V2 kr 

For r — > 0, obviously i p -*■ oo. This is due to the “turn-over” of the polar angle at r = 0 (compare 
the remark at the end of the problem). The e, -component of the acceleration is 



a r = r ■ e, = r — rtp 2 — —k I cos tp + sin <p ip I . 



2 kr 



2k 



The angular acceleration ip follows from 10.32, whereby v = 0: 

vr v 2 sin (p 

V ~ 2ry/lki ~ 4 r 2 

Equation 10.34 inserted in 10.33 yields 

D 2 ( 2 . 2 

a r — —k — - I r- cos (p + sin - tp 
4 r- \ k 

v 2 , 

= T (1 + 2 

4k - (1 + cos i p)~ ' 

3v 2 

a r = , radial acceleration. 

4 k 



ir 

2k 



cos <p + cos 2 1 p) — 



2k 



( 10 . 33 ) 



( 10 . 34 ) 



( 10 . 35 ) 



For the second component of the acceleration (azimuthal acceleration). 



a v = r • e p = rip + 2r<p 



v 2 sin i p 
4 r 



- 2k 



v 2 sinip 
2 kr 



3 v 2 sin i p 

4 r 



3 v 2 simp 

4 k 1 + cos i p 



( 10 . 36 ) 



Obviously, a v — > — oo for ip — > 180° (the angle <p turns over — compare the remark at the end of the 
problem). 

Because the acceleration components a r e r and a v t v are orthogonal, the magnitude of the acceler- 
ation is given by 



a = J# T^-= 3 - 1 1+ Sm2y = 3 - I 2 . 

V r v 4ky (1+ cos i p) 2 4 k y 1 + cos (p 

For the total acceleration, it also holds that a — > oo for (p -> 180°. 
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Remark: The angular velocity 

10.32 and the angular acceleration 
10.34 become infinite for r = 0. 
This singularity is implied by the 
choice of the coordinate frame and 
is independent of the motion along 
the cardioid. Consider, for exam- 
ple, the uniform motion of a parti- 
cle on an arbitrary path in polar co- 
ordinates. Let the origin be located 
on a normal to the path. Because 
co = v/d, the angular velocity de- 
pends on the separation between 
origin and path: a>\ < oh < . . . < 
co n . In the limit with the origin lo- 
cated on the path, the angular ve- 
locity becomes infinite. 




If the origin of the polar coordinates is located on the path, the 
angular velocity becomes infinite. 



Problem 10.2: Representation of a vector in cylindrical coordinates 

Write the vector A = zei + 2.re2 + ve 3 in cylindrical coordinates. 

Solution For the solution, we make the ansatz A = A e e e + A v e v + A,e,. The unit vectors of the Cartesian 
frame have to be replaced by those of the cylinder system. Moreover, the components, namely, z,2x, 
and y, have to be expressed by cylindrical coordinates. 

The system of equations 

e e = ei cos <p + e 2 sin <p, 
e 9 = — ei sin <p + e 2 cos <p 

may be solved for ei , and yields 

ei = e e cos <p — e,, simp, 
e 2 = e e simp + e v cosip. 

It further holds that 

x = gcoscp, y = psimp, z = z. 

Insertion yields 

A = z(e e cos <p — simp) + 2 q cosi p(e e siny) + cos <p) + q sin <pe z . 

Thus the components are 

A g = z cos cp + 2q cos i p simp, 

A v = 2 q cos 2 (p — z sin <p, 

A z = Q suup. 
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Problem 10.3: Angular velocity and radial acceleration 

A rod rotates about P\ in a plane with the angular velocity 
a> = ke smip . At the time t — 0, let cp = 0. The straight line 
intersects a fixed circle of radius a at the point P 2 . 

(a) Find the angular acceleration of the rod. 

(b) Find the velocity v e and the acceleration b e of the point 
Pi along the rod. 

(c) Find the velocity and the acceleration of the point P 2 with 
respect to the center of the circle. 

Solution (a) The angular velocity is 
to = <p = ke An<p 
=>■ for the angular acceleration 

d> = ip = fca>cos<pe sm? ’ 

= k 2 e 2sm,p cos cp . 

(b) The position vector to the point Pi on the rod is 
r = ge e , where q = 2a cos cp =P r = 2 a cos 1 p e e . 

The velocity of P 2 is obtained from the relation 

r = e e Q + . 

and the acceleration is 

f = (e ~ B<P 2 )e e + (QV + 2 W) e f> ■ 

Insertion yields 

v = f = —2a sin tp cp e e + 2a cos <p<pe v 
= 2 a (— sin tp <p e e + cos (p<pe v ) , 

b = r = 2a [(—95 sim/j — 2i p 2 cos cp) e e + ( ip cos< p — 2i p 1 sin^>) e,,] . 

For the velocity and acceleration along the rod axis, that is, in the e e -direction, one obtains 
v e = —2 a<p sin^> e e , b e = —2a (ip sini p + 2 cp 2 cos 1 p) e e . 

The negative sign indicates that both v e as well as b e point toward the center of rotation P l . 
(c ) The rotation angle of OP 2 equals 2<p, and the velocity along the circle is 

r Pl = a ■ 2cp => Vp 2 = 2acp = 2ake sm<p . 





Motion around P-\ with angular ve- 
locity u>. 



The normal acceleration is 
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the tangential acceleration is 



bp^ r = — — = 2 aip — 2ak 2 cos cp e 2sm<p , 
dt 

and the total acceleration is 






Vector Differential 
Operations 



Scalar fields: The notion of scalar field means a function <f>(x, y, z) that assigns a scalar, 

the value <p(x yi, zi), to any space point P(x i, y\, Zi). Examples are temperature fields 
T (x, y, z ) and density fields g(x, y, z) (e.g., mass density, charge density). 

Vector fields: A vector field correspondingly means a function A(.r, y, z) that assigns a 

vector A(jci , >'i , zi) to any space point P(x \ , y\, zi). 

Vector fields are, for instance, electric and magnetic fields, characterized by the field 
strength vectors E and H, or velocity fields v(x, y, z) in flowing liquids or gases. 



The operations gradient, divergence, and curl (rotation) 

Gradient: Given a scalar field <p(x, y, z) , the gradient of the scalar field at a fixed position 

Pq(xo, yo; ZoX denoted by grad <p(x o, yo, ZoX is a vector pointing along the steepest ascent 
of (p, the magnitude of which equals the change of </> per unit length of the path along the 
maximum ascent at the point Pq(xo, yo, ZoX 

In this way, any point of a scalar field can be associated with a gradient vector. The set 
of gradient vectors forms a vector field associated to the scalar field. Mathematically the 
so-defined vector field is given by the relation 

9 9 9 

A(x, y, z) = grad(/> = e!— <p + e 2 — 0 + e 3 — 0. (11.1) 

ox dy oz 

To simplify the mathematical description, the following notation is used: 

3 9 9 

grad 0 = V0, where V = d— + e 2 — + e 3 — • 

ox oy 9z 

(V : spoken “nabla” or “nabla operator”.) 
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VECTOR DIFFERENTIAL OPERATIONS 1 1 



Definition of an operator: The nabla operator is a symbolic vector (vector operator) 

that, when applied to a function p, generates the gradient of p. Taken as such, the operator 
is meaningless; it has to operate on something, for example a scalar function p(x, y. z ). 

We now demonstrate that the vector field V</> has the properties quoted above. For this 
purpose we need the total differential of p, namely 



3 p 3 p dtp 
dp — —dx + —dy + —dz. 
dx dy dz 



( 11 . 2 ) 



This quantity describes the main part of the total increase of the function p if x changes by 
dx, y by dy, z, by dz, that is. 



Ap (p{x + dx, y + dy, z + dz) — (fix, y, z). 



The Taylor expansion up to the first-order term yields 



0 (r + dr) — p(x + dx, y + dy, z + dz) 



d </> 30 30 

= p(x, y, z) + — dx + —dy + z—dz + ■ 
dx dy dz 



and therefore 



At/) = 0(r + dr) 



30 30 30 

0(r) = z~dx + —dy + —dz ■ 
dx dy dz 



— dip + terms of higher order. 



(11.3) 



This explains the name total differential for the main part of the total increase of the function 
t/>. We thereby have used the Taylor expansion of a function (up to the first terms in the 
small quantities dx, dy, dz). In Section 22 Taylor expansions will be outlined in detail and 
explained by numerous examples. We recommend that you have a look at this section now. 

Using the infinitesimal position vector dr — (dx, dy, dz), we may also write the total 
differential as follows: 



d(p — V0 • dr 



I 30 30 30 \ 
\ dx ’ dy ’ dz J 



(dx, dy, dz) 



90 , 9(/> 30 

= -z~dx + — dy + — dz . 
ox dy dz 



(11.4) 



Equipotential surfaces are surfaces on which the function cp takes a constant value, 
4>(x, y, z) = constant. 

As has been shown above, there is the relation 



V(p-dr — d(p, with dr — (dx, dy, dz). (1U5) 

Because dp represents the sum of the increases of p in each direction dr, dp = 0 means 
to stay on an equipotential surface. For this case, it holds that 



0 = dp — V(/> • dr E s, 



( 11 . 6 ) 
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where drEs lies in the equipotential surface 
ES. The scalar product V0 • dr es vanishes 
only then if the cosine of the enclosed angle 
vanishes (compare the opposite figure), pro- 
vided that V0 / 0. This implies that Vr/; and 
drEs are perpendicular to each other. Thus the 
gradient of (p is always perpendicular to the 
equipotential areas. 

We now consider the increase d(p along the 
gradient vector at a fixed point of the scalar 
field: Here dr is parallel to V(/> and then V r/j ■ d r 
takes the maximum value. Therefore, the vec- 
tor grad (p = V</> always points in the direction 
of the strongest increase of </>; see the opposite 
figure. 




Equipotential lines and the direction of the 
gradient. 



Divergence: Contrary to the gradient operation, the divergence is applied to vector fields. 

Given a vector field A = (A x , A y ,A z ), we further imagine a cuboid-shaped “control 
volume” (rectangular box) with the edge lengths At, Ay, A". 

The “vector flow” across an area represents the entity of vectors penetrating it perpen- 
dicularly, that is, the normal components of the vectors integrated over the entire area. 

The lateral faces of the cuboid are denoted by 

*^1 5 • • • 5 $ 6 - 

We now calculate the vector flow across all 
lateral faces of the cuboid (rectangular box). 

The edge lengths At, Ay, A" shall be chosen 
so small that the vector on the cuboid faces may 
be considered as nearly constant, such that the 
integration of the vector across the faces may be 
replaced by a simple summation. We shall count 
the vector flow as positive if it flows out of the 
volume, and negative if it flows into the volume. 

The vector flow through the faces is 

v i : -A a (t)AvAz, 

S 2 : A x (t + At) Ay At , 

53: — A y (y)ATAz, 

.sq: A y {y + Ay)ATAz, 

and in the third space direction is 



1 k 

i 

y+Ay 


‘A y (y+Ay)e } , 


Aflx)e x 

► 

y 


5 54 

5 3 


AJx+Ax)e x 

► 

y 

► 


a 


i 


A v (y)e 



t+At 



The flow across a cuboid. The extension in 
the z-direction (out of the paper plane) is 
not shown. 



S 5 : — A z (z)Ax Ay, 

s 6 : A z (z + At) At Ay. 



( 11 . 8 ) 
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A Taylor series expansion up to terms of first order, which is satisfied for small Ax, Ay, A z, 
yields 

9 

AAx + Ax, v. z) = A , (x . v, z) H AAx. v. z)Ax H , 



A z (x , y, z + A z) — A z (x, y, z) + — A z (x, y, z)Az H . 

9z 

The terms indicated by dots • • ■ are of higher order in the small increments Ax, Ay, Az 
and may be neglected. The resulting vector flow through the control volume follows by 
summation over the lateral faces: 

(A*(x + Ax, y, z) - A x (x, y, z)) AyAz 

+ (A y (x, y + Ay, z) - A y (x, y, z)) AxAz 

+ (A z (x, y, z + Az) — A z (x, y, z)) Ay Ax, 

9 9 9 

= —A x (x, y, z)AxAyAz + — A y (x, y, z) Ax AyAz + — A z (x, y, z)AxAyAz 



Thus the “flow” (total flow) through an infinitesimally small volume (Ax —>■ dx. Ay — » dy, 
Az dz) reads 



This last relation may be interpreted as analytic definition. As as flow of the vector field 




( 11 . 9 ) 





( 11 . 10 ) 



The expression in brackets is called divergence of the vector 
field A: 





( 11 . 11 ) 



Thus, the divergence represents the vector flow through a 
volume AT per unit volume. It may also be written in the 
form 



div A = V ■ A(x, y, z). 



( 11 . 12 ) 



Illustration of the divergence 



has been shown, it is identical with the geometric definition, through a volume, 
namely: 



div A = lim 

A V— >0 



flow of the vector field A through A V 



AV 




While the argument of the gradient operation is a scalar, the divergence represents the 
scalar product of the operator V and the vector A. For a vanishing divergence, the total flow 
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through an infinitesimal volume equals zero, that is, the in-flow just balances the out-flow. 
If at some point of the vector field div A > 0, one says that the vector field there has a 
source; for div A < 0, one speaks of a sink of the vector field. This is immediately clear 
from the definition of the divergence as net flow = out-flow — in-flow per unit volume. 

Curl (Rotation): 1 The operation curl A assigns a vector field curl A to a given vector 
field A. The vector field curl A informs about possible “vortices” of the field A (a vortex 
exists if there is a closed curve in the vector field fulfilling the condition that the contour 
integral § A • ds ^ 0 — see theorem of Stokes). The mathematical formulation of curl A is 
given by 

1 . curl A = V x A, or 

2. n • curl A = lim A/ .^ 0 A • d s) /A F. 
n is a unit normal vector on A F. 

The second definition states that the rotation may also be determined 
by forming the contour integral. The integration is performed over the 
vector field along a curve. More strictly speaking: One integrates over 
the projection of A onto ds along the tangent to the curve forming the 
border of A F. After division by A F, this yields the component of curl A 
along n. 

The rotation is thus determined by two distinct definitions. The first 
detail 




Oriented surface 
element. 

of these reads in 



ei 



curl A = V x A = 



9/9x 

A x 



e 2 e 3 

d/dy a/9z 
A,, A- 



= ei 







+ e 3 





( 11 . 14 ) 



One has to prove that both definitions are iden- 
tical. Here we show the identity only for the 
x-component. 

x -component of the curl of A: One may 

integrate about an area A F = 4AyAz in the 
y, z-plane (see the lower figure), n then points 
along the x-axis, that is, n • curl A just yields 
the x-component of curl A, namely (curl A)*. 




Illustration of a vector field A with vorticity on 
surface element A F with normal vector n. 



1 In German literature, curl is named rotation (rot), so curl A = rot A. 
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For the loop integral one has 




c c 



f 



+ A-dz, 



f 



— fh(A y dy + A z dz), 



(11.15) 



since for this orientation of the area concerned (see figure) Ax = 0 (i.e.. At does not enter 
at all). In other words: Because x remains unchanged (dx = 0), f A x dx drops out. For the 
exact definition of the contour or loop integral we refer to Chapter 12. It is recommended 
to study this section in brief right now. 



Remark: § shall indicate that the integration is performed over a closed curve (contour 

or loop integral) in the counter-clockwise direction. j : - means integration over a section of 
the curve. For calculating a contour integral, we employ the values of the functions in the 
middle of the individual sections (marked points). 




Pi Pi A Pi 

hhhf 



( Aydy + A z dz) 



& A y (x, y,z — Az)2A_y + A.(x, y + Ay, z)2Az 
— A y (x, y, z+ Az)2Ay - A z (x, y — Ay, z)2Az. 
According to the Taylor expansion this yields 



(11.16) 



dAy 1 

A > - nt Az 



2 Ay ■ 



dA- 1 
— ^Ay 

dy J 



2A z 
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9 Ay 

3 z 



■A z 



2 Ay — 



dA. 

A - -^~ A y 

3 y 



2A z 



= AAyAz 



dA 

3 y 



8 Ay 

3 z 



The enclosed area is AT = 4 Ay At. From that, it follows that 



lim 

Af^O 



A ■ ds\ 3 A z 3 A v 

A F ) x dy dz ' 



(11.17) 



Hence, the x -components corresponding to the two definitions of curl A coincide. For 
the remaining two components, the equivalence of the definitions may be demonstrated in 
an analogous way (which will be skipped), q.e.d. 

The second definition of curl A implies that the rotation at some point of the field A 
vanishes if the contour integral (f A ■ ds (loop integral) enclosing this point equals zero — 
see the theorem of Stokes. From there originates the name “rotation.” A finite value of the 
loop integral expresses a certain rotation, that is, vortex formation of the vector field (to be 
visualized as a flow field). 



Multiple application of the vector operator nabla: Given a scalar field /( r) and a 

vector field g(r), then 
(a) 

3 2 f 3 2 f 3 2 f 

V • (V/) = —A + —A + — y = divgrad f(x,y,z) = A f(x, y,z), (11.18) 

ox z dy z az 

where A is introduced as a new operator: 

3 2 3 2 3 2 

A = — -H -H = V • V 

dx 2 3y 2 3 z 2 



(A, spoken: delta, is called the Laplace operator 2 ). V • (V /) = div V/ is a scalar field. 



2 Pierre Simon Laplace , b. March 23, 1749, Beaumont-en-Auge — d. March 5, 1827, Paris. After his school 
education Laplace became a teacher in Beaumont and, by mediation of D’Alembert, became appointed as professor 
at the Military School of Paris. Because Laplace used to quickly modify his political convictions, he was swamped 
with honors both by Napoleon and by Louis XVIII. Among his works his Analytic Theory of Probability (1812) 
and the Celestial Mechanics (1799-1825) became significant. The theory of probability calculus contains, for 
example, the method of the generating functions, the Laplace transformations and the final formulation of the 
mechanical materialism. The Celestial Mechanics presents, for instance, the cosmologic hypothesis of Laplace, 
the theories of the earth’s shape and of the moon’s motion, the perturbation theory of planets, and the potential 
theory with the Laplace equation. 
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(b) 



V x (V /) = curl grad / = 



ei 


e 2 


e 3 


_a_ 


d_ 


3 


dx 


dy 


Bz 


§L 


df 


a if 


dx 


dy 


Bz 



= 0. 



Thereby it is of course required that / is twofold continuously differentiable. The physicist 
always presupposes functions that are sufficiently often continuously differentiable; this is 
also assumed below. Hence, a gradient field has no vortices! 

(c) 



V (V • g) = V 



V dx dy dz ) 



9 ( dgx dgy dg : 

dx \ dx dy dz 



, 9 / 9gx , 3 gy , 3 gz 

ei + 9V ( 37 + + ^7 » 62 



3 / dg x dg y 

dz \ dx dy 

= grad(div g) is a vector field. 



3 gz 
dz 



e 3 



(d) 



V • (V x g) = divfcurl g) = 0. 

Hence, a rotation field has neither sources nor sinks, as 
is graphically clear: The vector field A = co x r with 

co =constant is so to speak an optimum vortex field (the 
velocity field of a rigid body rotating with the angular 
velocity co). 

For this maximum vortex field one has curl A = 2co. that is, it is a 
constant vector field that obviously is divergence-free. One should 
note the similarity of V with a vector: The triple scalar product 
involving identical vectors vanishes. 

(e) 

V x (V x g) = curl(curl g) 

= V(V • g) - (V ■ V)g 
= grad(div g) - Ag 




to 



The velocity field of a rotating 
rigid body: A = »xr. 



The rotation of the ve- 
locity field above: curl A = 
curl co x r = 2co. 



is a vector field. 

The proof is simple, because according to the expansion theorem, 
C x (B x A) = B(C • A) - (C • B)A. 



This twofold application of the rotation operator physically and geometrically means that 
the vortices of the vortex field are calculated. 
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(f) 

div(B x C) = C • (curlB) - B • (curl C). 

(g) 

v-(/g) = V/.g + /V-g. 

(h) 

v X (/g) = V/ x g + / V x g. 



Problem 1 1 .1 : Gradient of a scalar field 

Given the scalar field ip = x 2 + y 2 = r 2 , find the gradient of ip. 

Solution 

Vi p = 2(xe v + yt y ) = 2 sj x 2 + y 2 e r = 2re r . 



Problem 1 1 .2: Determination of the scalar field from the associated gradient field 

Let V<p = (1 + 2xy) e v + (x 2 + 3 y 2 ) e y . Find the associated scalar field. 



Solution 



d <p 
dx 
dip 
dy 



= (1 + 2xy) 
= (x 2 + 3y 2 ) 



By comparison: 
fi(y) = y 3 + c 1 . 



=> v>(x. y) = x + x 2 y + fi(y), 

=J> <p(x, y) = x 2 y + y 3 + f 2 (x). 



f 2 (x) = jc + C 2 ; 



thus, 

ip{x, y) = x + x 2 y + y 3 + C . 



Problem 1 1 .3: Divergence of a vector field 

Calculate the divergence of the field of the po- 
sition vectors: 



Solution 



r = xei + y e 2 +te 3 . 



dx dy dz , 

divr = 1 — - -\ = 3. 

3x dy dz 



Thus, the vector field r everywhere has a 
finite divergence (i.e., source density) of mag- 
nitude 3. To generate this field in practice by a 
flow, one would have to attach sources of inten- 
sity 3 to any space point. 




The field of the position vector A(x, y, z) = r. Flow 
out of the volume element A V is larger than flow 
into the volume element. 
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Problem 1 1 .4: Rotation of a vector field 

Calculate the rotation of the vector field 



A = 3x 2 yei + yz 2 e 2 - xze 3 . 



Solution 



curl A = ei 



3 (-xz) 3 yz 2 

dy 3 z 



d3x 2 y 3 (—xz) 



3 yz 2 d3x 2 y 

dx dy 



= —2 yz ei + z e 2 - 3x 2 e 3 . 



Problem 1 1 .5: Electric field strength, electric potential 

Let a positive electric charge of magnitude Q be localized at the origin of the coordinate frame. The 
field intensity E describing the electrostatic field is given by 

E= ^e r . 



Solution 



where r denotes the spatial distance front the coordinate origin, and e, represents the corresponding 
unit vector in radial direction. Calculate the associated potential field (let U denote the potential field, 
then E = —VIZ) and show that it satisfies the Laplace equation A U = 0, except for the origin. 



Q Qr 

E = — i-e, = -r- E = —V U. 



(11.19) 



Because E points in the radial direction and the gradient means the derivative along this direction, 
one has 

dU 

m = -—, 

dr 

and because E is a function of r only, it follows that 

( f dr 1 

U = ~ |E | dr = -Q / - = Q- + C. 



One easily confirms the relation 11.19 for this potential field, for example, for the x -component 



dU 3 Q 



dr d /I 



= — =~Q \ - ) = Q— = E 



dx dx r 



dx dr V r 



etc. The constant C is usually set to zero, that is, the potential vanishes for r — > oo. 

| 3 x 3 y 

{ dx ( x 2 + y 2 + z 2 ) 3/2 3y (x 2 + y 2 + z 2 ) 3/2 



( 11 . 20 ) 



3 z (x 2 + y 2 + z 2 ) 3/2 



( 11 . 21 ) 



div(V(7) = AU — 0 for r ^ 0. At r = 0, one has divE = — div VI/ ^ 0 (see below: the Gauss 
theorem). 
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Problem 11.6: Differential operations in spherical coordinates 

Given a scalar field 0(r, ft, ip) and a vector field A (r, ft, <p), 
which are the relations for (a) V0, (b) V • A, (c) V x A, 

(d) V 2 0 in spherical coordinates? 



Solution (a) Gradient: For the total differential, it holds that 



d(f> = V0 ■ dr. 

In spherical coordinates one has 

3 0 30 3 <b 

d< b = -rr-dr + -^-dft + -a - dtp , 
or oft dip 

dr dr dr 

dr = — dr H dft H dtp 

3 r dft dep 



= e r dr + re# dft + r sin ft e v dip. 



( 11 . 22 ) 



(11.23) 




(1 1 .24) Spherical coordinates: A reminder. 



and 



V0 = (V0),.e,. + (V0)#e# + (' V<b)<p£<p ■ 

The partial derivatives of the position vector have already been calculated in Chapter 10 when 
formulating the unit vectors: 



3r 



3r 



— = e r , — = re#, 
dr dft 



dr . 

— = r sm ft e„ . 
dip 



(11.25) 



By insertion and comparison of coefficients, 11.22 immediately yields for the components of the 
gradient in spherical coordinates 



30 30 30 

CVtp) r dr + (V0)#r dft + (V0)„r sin ft dip = — dr H • dft -\ dip, 

dr dft dip 



30 1 30 1 30 

V 0 — — e,. + e# T Cm 

dr r dft r sin ft dip 

= (V0) r e r + (V0)#e# + (V0) ¥ ,e v . 



(11.26) 



(b) Divergence: The divergence may be expressed by the flow of the vector A across the surface 
of an infinitesimal volume element A V : 



divA = 



lim 

A V-»0 



/af A • ndF 
AV 



(11.27) 



The figure shows the volume element with the magnitude 

AV = r 2 sin# Ar Aft Atp. (11.28) 

Calculation of the flow components (to first approximation): The flow in the e, -direction across 
the area ADHE is 



A(r, ft, ip)e r AF r = A, r 2 sin ft Atp Aft. 
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the flow across the back area BCG F is 

A(r + A r, #, <p)e, r AF r+Ar = A r r 2 sin # A <p A# + — (r 2 sin # A,. Aip A;?) A r. 

dr V ' 

The difference yields the contribution of the flow in the e, -direction to the surface integral in 1 1.27. 
The flow excess is 

sin# — (r 2 A r )A<p A# Ar. (11.29) 

dr 

The flow excess in the e# -direction (areas ABFE and DCGH) correspondingly follows as 



r — (sintfA,)) A <p Ar A#. 
3# 



(11.30) 



The flow excess in the e„ -direction is 



3 

r — A v Ar A# A<p. 
dcp 



(11.31) 



Summation of the contributions 11.29, 11.30, 11.31 yields the flow integral <f A • n dF. Then, 11.27 
yields the expression for the divergence: 



13 , 

V • A = — — {r 2 A r ) — 



1 



-(sin# A#) + 



1 3 

Ain . 

r sin # 3i p 



(11.32) 



r 2 dr ' r sin # 3# ' 

(c) Curl (Rotation): The geometric definition traces the rotation operation back to a contour integral: 
A ■ ds 



n ■ curl A = lim 



AF — *-0 A F 



(11.33) 



Component along e r : 

The e r -component of the rotation is obtained when performing the contour integral along the curve 
ADH EA( n = e, ). The enclosed area is then 



AF = r~ sin # A# A<p (compare fig.), 

D H E A 



(11.34) 



ADHEA 



A • ds — J* + T j* T . 



A D H E 

The partial integrals are 

D 

A • ds = A • e#r A# = A^r A#, 



/- 

A 

A 

/ 



A • ds = A • (— e v )r sin # Ac -p = —A v r sin # A <p. 
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And to first approximation, 

E 



A • ds = — ^rA^A)? + — (rA#A$)AipJ , 



A • ds = r sin i?A v A<p H (r sin dA v Aip)Ati. 

d& 



Then the contour integral along the closed curve is 

/ 9 9 

< p A • ds = r — (sin;? A v )Aip A) ? — r — (A#)A& A<p. 

Jadhea 30 9 <p 

From 1 1.33, 1 1.34, 11.35, it follows that the e r -component of the rotation is 

1 r 9 a . 1 



curl A = 



1 

r sin ]? ]_i 



-(sint?A,,) - —A# 
dip 



Accordingly for the curve AEFBA with A F = r sin d Ar A<p, 

/a a 

( p A ■ ds = ( A v r smtiA<p)Ar H (A, Ar)Atp 

Jaefba 9 r 9 ip 

and because n = e#, it follows that 



curl A = 

i» r sin 



— — I" — A r 

mt? 9 < p 



- sin)? — (rAp) 
dr 



Investigation of the curve ABC DA yields 



curl A = 



(r Aff) — ——A, 



<p r \dr did ) 

The results 11.36, 11.37, 1 1.38 may be combined into a determinant: 



V x A = 



(d) Laplace operator: The Laplace operator is defined by 
V 2 ^> = div Vtp. 

Using the results 11.26 and 11.32, it follows that 



e r 


re & 


r sin ft 


a 


d 


_d_ 


dr 


9# 


dtp 


A r 


rA & 


r sin )? A v 



1 3 / , 3 <p 



V 2 </> = — r 



r 2 dr \ dr ) r 2 sin & dd \ 9 1 ? / r 2 sin 2 ]? dip 2 



i a 2 ^ 



(11.35) 

(11.36) 

(11.37) 

(11.38) 

(11.39) 

(11.40) 
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Differential operators in arbitrary general (curvilinear) coordinates 

In Chapter 10 we outlined curvilinear coordinates (e.g., spherical and cylindrical coordi- 
nates). In Problem 11.6, the differential operators V, div, and curl have been derived in 
spherical coordinates, basing on special considerations. Now we shall develop the general 
approaches for calculating differential operators in arbitrary curvilinear coordinates. 



Brief repetition: Let r(x v ) = x v e„ be the position vector in Cartesian coordinates 

x v (v = 1,2, 3) that are related to the curvilinear coordinates q a (a = 1, 2, 3) via x v = x v 
(c/i , qi, t/ 3 ) . The x v may then be inserted in the position vector, which yields 



r(x v ) = r (x v (q a )) = r (q a ). 



(11.41) 



The new unit vectors e ?[r , which in general are characteristic for the point q„. may be 
defined at each point q a (er — 1, 2, 3): 



9r {q lx )/dq a 

|3r(^)/a^| ’ 



cr = 1,2,3, 



3r (gy) 

3^cr 






with h a 



3r (g M ) 

dq<r 



(11.42) 



(11.43) 



Here the h a (a = 1,2, 3) are scaling factors. The unit vectors e 9(j point along the q a - 
coordinate line toward increasing q„ . 

The coordinate areas are obtained by solving the three equations x v = x v (q l , t/ 2 , l/s) 
for q a : 



q a = q n (x\,x 2 , x 3 ) = q a (x IJL ) . 



(11.44) 



q a = constant = c n - (o — 1 , 2, 3) are the equations for the coordinate areas. 

One may now construct other unit vectors 
E qo at the point P(x, y, z) — P(qi , qi) 

(see the figure), namely 

Vq« 



— 



|V<?a 



a = 1,2,3. 



(11.45) 




U, 



Different basis vectors at the same point P. 



The E Ch are obviously perpendicular to 
the coordinate areas q a = c„- . Thus, there 
are two sets of unit vectors at each point 
P (q a ), namely e ?tr and \i, h . In general, these 
sets are distinct. We shall demonstrate in the following that these two basic frames coincide 
only then if the curvilinear coordinates are orthogonal. One also has to take into account 
that both e q „(q\, qi, qi) as well as E ?(j (gi, qi, qi) depend on the point P(q\, qn, qi), that 
is, their orientations in general vary from point to point. 
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An arbitrary vector A may now be expressed in terms of both the base e g<j as well as 
the E qa , 

A = Aie 9l + A 2 e 92 + A 3 e ?3 = (?iE gi + a 2 E q2 + a 3 K q3 . ( 11 . 46 ) 



The A, and a,-, respectively, are the components of A in the bases concerned. Instead of the 
normalized base vectors e, la or E, /n the nonnormalized vectors 



K 

and 



dr(q a ) 

3 q v 



(v= 1,2,3) 



B g> = V 9v (v= 1,2,3) 



( 11 . 47 ) 

( 11 . 48 ) 



may also be used. They are called unitary base vectors and are in general not unit vectors. 
For an arbitrary vector A, 

A = Ci- b C 2 — b C 3 - — = Cib 9l + C 2 b 92 + Cib q3 ( 11 . 49 ) 

oq\ dq 2 dq 3 

and 

A = C\Wq\ + c 2 V q 2 + c 2 Vq 2 = ciB gi + c 2 B g , + c 3 B g3 . 



The components C v ( v = 1, 2, 3) are called contravariant components and c v (v — 1, 2, 3) 
covariant components of the vector A. They play an important role in the general theory of 
relativity where all coordinate frames are used on equal footing. In Cartesian coordinates 
the co- and contravariant components of a vector are equal to each other, as is immediately 
clear from their construction. 

(a) Arc length and volume element: From r = r(q \ , q 2 , q 2 ), one obtains 



3r 3r 

dr — - — dq 1 + - — dq 2 + 
oqi dq 2 



3r 

—dq 2 = hidqi e 9l + h 2 clq 2 e g , + h 3 dq 3 e 93 
oq 3 



( 11 . 50 ) 



Therefore, for the differential ds of the arc length, it results that 

ds 2 — dr ■ dr, ( 11 . 51 ) 



which for orthogonal coordinates (e q/J ■ e gv = S^ lv ) simplifies to 

ds 2 — h 2 x dq 2 + h\ dq\ + h\dq 3 . ( 11 . 52 ) 

For nonorthogonal coordinates, it holds that 

b„. ; • b qv — 1 1 u h • e g „ = gjiv 7 ^ h jj h yl) // v , (11.53) 

and therefore it follows from (1 1.50), (1 1.51), and ( 1 1.53) that 
(ds) 2 = dr ■ dr 

\ 

v- , 

— T, g t ivdqi,dq v . 

H,v 





( 11 . 54 ) 
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This is the fundamental quadratic (or met- 
ric) form. The g^ lv are called metric coeffi- 
cients (since they determine the measure- 
ment in the coordinates q v via the length 
element ds 2 ) or also metric tensor (briefly: 
metric). If — 0 for /i ^ v. the coordinate 
frame is orthogonal. In this case gn = h\, 
g 2 2 = h\, ^33 = h\. The metric tensor is 
of basic importance in the general theory 
of relativity. It is determined there from the 
energy (mass-) distribution in space. 

The equations enabling this are called f /i 

Einstein equations. Illustration of the volume element. 

The volume element dV may easily be calculated for orthogonal coordinates (see the 
figure): 

dV = | (/? i dqi e 9l ) • [(h 2 dq 2 e qi ) x (/t 3 r%e ?3 )] \ 

= h \h 2 h 2 dqi dq 2 dq 2 , (11.55) 

because 

| e 9 l ' ( e <?2 X e <? 3 ) | = 1 - 




Problem 1 1 .7: Reciprocal trihedral 

Given are the three noncoplanar vectors a. b, c for which a-(bxc) 7 ^ 0. Show that the three reciprocal 
vectors 



b x c 

a ■ (b x c) ’ 



b' = 



a • (b x c) 



a x b 
a ■ (b x c) 



(11.56) 



are also noncoplanar and that 
(a) 



a' • a = b' ■ b = c' • c = 1 , 

(b) 

a' • b = a' • c = 0, 
b' • a = b' • c = 0, 

c' • a = c' ■ b = 0. (11.57) 

(c) If a • (b x c) = V, then a' ■ (b' x c') = 1/ V . 



a • (b x c) 

a = = 1. 

a • (b x c) 



Solution (a) 
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In the same way one may conclude that 
b' • b = c' • c = 1 . 



(b) 

, b(bxc) 

a • b = b • a = = 0, 

a • (b x c) 

and similarly for the other cases. 

(c) One has 

, b x c , c x a 
b = 



a x b 

V 



V V 

Then it follows that 

(b x c) • [(c x a) x (a x b)] 

(a x b) • [(b x c) x (c x a)] 

V 3 



a' • (b' x c') = 



(a x b) • [c • ((b x c) • a) — a • ((b x c) • c)] 

F 3 

[(a x b) • c][(b x c) • a] [a • (b x c)] 2 V 2 
V 3 = F 3 = F 3 



1 

V' 



From there it follows that a', b', c' are noncoplanar if a. b, c are noncoplanar. 



Problem 1 1 .8: Reciprocal coordinate frames 

Let qi, q 2 , q 3 be general coordinates. Show that dr/dq\ , 3r/3 q 2 , 3r/3 q 3 and V^, Vq 2 , Vq 3 form 
two reciprocal systems of vectors and that 

3r / 3r 3r \ 1 

7- • hp- x r [Vqi ' (V?2 x v?3)) = 1 • 

oq l \oq 2 9^3/ J 



Solution One has to show that 

3r 1 1 for v = ti, 

‘ — l 

[0 for v / /i , 

where /j., v may take any of the values 1,2,3- Now 

3r 3r Sl- 

ur — - — dq\ + - — dq 2 + - — dq 2 . 
oq i oq 2 oq 3 

and therefore after multiplication by Vq l 

( 3r \ ( 3r \ ( 3r \ 

Vqi ■ dr = dqi = I V^i • — J dq\ + I Vq { • — J dq 2 + I • — J dq 3 . 



( 11 . 58 ) 



From there it follows that 
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The other relations result in a similar way by forming Vq 2 ■ dr — dq 2 and V</ 3 ■ dr = dq 2 . Thus, 
the reciprocity of the vector systems dr/dq v and Vq v is demonstrated. 

From the preceding problem it then immediately follows that 

3r / 3r 3r \ 1 

T " ( x — ) 1 {V 9l • (Vq 2 x V<? 3 )} = 1. 

dqi \dq 2 dq 3 ) \ 

This statement is equivalent to the following theorem on Jacobi determinants: 



V X,y,z ) 



def. 



V<y, • (V q 2 x Vg 3 ) = 



dqi 


dqi 


dqi 


dx 


dy 


dz 


dq 2 


dq 2 


dqi 


dx 


dy 


dz 


dq 3 


dqi 


dqi 


dx 


dy 


dz 



which reads 



j( X ’ y ' Z )-j( qi - q2 ’ q3 ) = l. 

V?1^2,?3/ V x,y,z ) 

One can check easily that 3r/3<j>, and V</, fulfill relation 1 1.56 from Problem 11.7. 



(11.59) 



(b) Gradient in general orthogonal coordinates: Let tp(q \ , q 2 , qz) be an arbitrary 

function. We look for the components f \ , f 2 , f 3 of the gradient in the general base e cjv , 
that is. 



= fltqi + f2^q 2 + h^q 3 - 

Because 

3r 3r 3r 

dr = — — dcji + — dq 2 + —dq 3 
dqi dq 2 dq 3 

= h\e qi dqi + h 2 e qi dq 2 + h 3 t m dq 3 , 
because of the presupposed orthogonality of the e ?v , it follows that 
dtp — Vtp ■ dr — hi f i dqi + h 2 f 2 dq 2 + 773/3 dq 2 . 



But it also holds that 

dtp dtp dtp 

dtp — —dqi + — — dq 2 + —dq 3 . 
dqi dq 2 dq 3 

A comparison of the last two relations yields 



Sj(p — ^ ^ ^ 

h\ dqi h 2 dq 2 h 3 dq 3 

In operator notation this reads 



V = 



1 



hi dqi ) + 692 (/? 2 3ft) + tq ' Ol 3 <73 ) ' 



(11.60) 



(11.61) 
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From there it follows especially for cf> — q\ that 

V?! = (11.62) 

hi 



and therefore |V^i| = l//zi or generally |Vz/„| = 1 /h v (v = 1, 2, 3). 
Because 



E* = 



Vgv 

|V? v | 



(compare to (1 1.45)), it results that 



E a = 

|V 9v 



Vo y 

= = h v — - — e 9v 

h, v 



(v = 1,2,3). 



This means that/or orthogonal coordinates the reciprocal base systems E ?v and e 9i , coincide. 
This happens, of course, in particular for Cartesian coordinates. 

For the following, the relations 



e 9l = h 2 h 3 Vq 2 x V^ 3 , 
e ?2 = h 3 h\Vq 3 x Vgi, 

e ?3 = hih^qi x V^ 2 (11.63) 

are helpful. They may be checked quickly, for example, 

h 2 hiVq 2 x V<? 3 = /z 2 /7 3 ( ^ x ^ ) = ^(e ?2 x e 93 ) = e„. (11.64) 

\«2 «3 ) h 2 h 3 - 

(c) Divergence in general orthogonal coordinates: We shall now calculate 

divA = V • (Aie 9l + A 2 e qi + A 3 e 93 ) 
in general coordinates. For this purpose we consider at first 



V • (Aie ? i) = V • ( Aih 2 hiWq 2 x V<? 3 ) 

= (V(Ai/z 2 /? 3 )) • (V q 2 x V< 7 3 ) + A x h 2 h 3 V • (V^ 2 x V<? 3 ) 



= (V(A 1 /z 2 /7 3 )) • ( ^x^l+0 

n 2 h 



= V(Ai/j 2 /7 3 ) 



/7 2 /73 



• j — t — {A\h 2 h 3 ) + -p- — (A l h. 2 h 3 ) + -p- — (A\h 2 h 3 ) 
L«t 9<?i « 2 9g 2 773 9^3 

1 9 

" (A 1 /7 2 /z 3 ) . 



h 2 h 3 



h\h 2 h 3 dq 
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Similarly it follows that 

1 9 

V • (A 2 e 92 ) = 
and 

V • (A 3 e 93 ) = 

Therefore, 



h\h 2 h 3 9<?2 

1 9 

h i /z 2 /t 3 9<73 



(A 2 hih 3 ) 



(A 3 h\h 2 ). 



divA = V • (Aie 9l + A 2 e q2 + A 3 e ?3 ) 

= V • A ] , + V • A 2 e 92 + V • A 3 e 93 , 

1 r 9 9 9 

dlvA = , , 7 - — (Ai/t 2 ^3) + — (A 2 h l h 3 ) + — (A 3 /z i/t 2 ) 

h\h2rl2 \_Otfl OCfa 9(^3 



(11.65) 



(11.66) 



(d) Curl (Rotation) in general orthogonal coordinates: We have to calculate 

V x A = V x (Aie ?1 + A 2 e 92 + A 3 e 93 ) 

= V x (Aie 9l ) + V x (A 2 e 92 ) + V x (A 3 e ?3 ) . 

It suffices to consider, for example, the term V x ( A 3 e^, ) in more detail. We obtain 

V x (Aie 9l ) = V x (Ai/?iV<?i) = V (A i /z i ) x V^i + A\h\ V x Vq\ 

= V(A 1 A 1 )x ^+0 
h i 

e„ 9 e fl2 9 e a3 9 

■7^7— (Ai/n) + — (Aj/m) + — (At/zt) 



L /z i 9?1 
h 3 hi dq 3 



h 2 dq 2 



h 3 dq 3 



(AM - 



9 



h\h 2 dq 2 



(AM. 



Therefore, 



V x A = 






, . | t — (A 3 /j 3 ) - - — (A 2 h 2 ) 
h 2 h 3 ldq 2 dq 3 






A (Ai/?i) _ JL (A3 / ?3 ) 

“3^1 |_9g 3 



9 9 

— (A 2 /t 2 ) - — (AiAi) 
oq i 9g 2 



h\h 2 

In determinant notation this reads 





h\t qi 




h 3 e q} 


1 


9 


9 ' 


9 


h\h 2 h 3 


9<?i 


dq 2 


9^3 




A\h\ 


A 2 h 2 


A 3 h 3 



(11.67) 
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(e) The delta (Laplace) operator in general (orthogonal) coordinates: One has to 

calculate Ai jr in orthogonal curvilinear coordinates. This does not provide any difficulties, 
because 

hi dq\ h 2 dq 2 h 3 dq 3 J 



A\ jr = V ■ Vi /f = V 
Using now equation (1 1.66) for the divergence, where obviously 



1 3 

Av=|-t- (v = 1, 2, 3), 
h v oQv 

one immediately finds that 
Ai J/ — V • Vi j/ 



( 11 . 68 ) 



1 



h\h 2 h 3 



3 / h 2 h 3 3 \[r 

dqi V hi 3 qi 



3 / h 3 hi 3i j/ 

dq 2 \ h 2 dq 2 



3 / h \h 2 d\// 
dq 3 V h 3 dq 3 



(f) Examples of special orthogonal coordinate frames 

1. Cylinder coordinates 

r(x, y, z ) = xei + yea + ze 3 

= pcos^oei +psini/)e2 + ,-e 3 = r (g,<p,z). 

Here q > 0, 0 < q> < 2n, — oo < z < oo. 

We identify q i — q, q 2 — <p, q 3 = z. According to (1 1.42), it then follows that 

e 9l = e e = cos cp ei + sin ip e 2 , 
e ?2 = e^, = — sin <p ei + cos cp e 3 , 

e ?3 — = ^3- 

Moreover, 
hi = h B = 



h 2 = /?„ = 



(11.69) 



(11.70) 



h 3 = h- 



3r 


= 1 , 


~dQ 


3r 




dip 


= e> 


3r 


= i . 


3 z 



According to equation (1 1.61), it therefore follows that 




